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ABSTRACT 


A Solar Power Satellite (SPS) microwave power density of 
20 mW/cm2 is confirmed as the level where nonlinear inter- 
actions may occur in the ionosphere, particularly at 100 
hm altitude: Radio wave heating at 100 1cm at this level 

can be produced in a test program outlined and the effects 
observed. A crucial element of the test program is a set 
of diagnostic measurements to assure that the frequency 
scaling used in the test is understood and that the results 
of the test can be extended to the SPS frequency. 

The limited test conducted in June 1977 at the Arecibo 
Observatory produced negative results for radio wave heat- 
ing of an underdense ionosphere, the condition to be en- 
countered by an SPS system. The results are consistent 
with available theory. Overdense heating produced stria- 
tions in the ionosphere easily detected nearly 700 1cm away 
by a small radar located to satisfy the predicted aspect 
sensitivity of the targets, as well as by an incoherent 
backscatter radar at Arecibo. These observations contain 
new results supporting some theories and casting doubt on 
others . 

An assessment of the radio frequency interference problem 
from striations shows that under some conditions it might 
be severe, but this can be established in the test program 
proposed. The effects of thermal self- focusing originally 
thought to be a potentially serious source of radio fre- 
quency interference are shown to be limited severely geo- 
graphically. The aspect sensitivity of field-aligned stria^ 
tions, regardless of how they are produced, makes interfer- 
ence-free regions above magnetic latitude about 60°. 

A test program is proposed to simulate the interaction of 
the SPS beam with the ionosphere, to measure the effects of 
the interaction on the ionosphere and on communication and 
navigation systems, and to interpret the results. 
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1. RESULTS 


1.1 NONLINEAR INTERACTION (See Section 3.1, 3.2) 

Theoretical studies confirm the prediction of Meltz 
that a 5 GW Solar Power Satellite (SPS) with power 
densities of 20 mW/cm? will heat the ionosphere's 
D, E and P regions sufficiently to excite nonlinear 
effects . The same and higher levels of radio wave 
heating can be introduced into the, ionosphere by high 
frequency (HP) waves and the effects observed (see 
Section 5 ) . 


1.2 LIMITED TEST PROGRAM (See Section 4) 

A limited test using the existing equipment at Arecibo 
at HP, 430 MHz and S-band for radio wave heating of the 
underdense ionosphere produced no striations over Are- 
cibo that could be detected with a sensitive radar at 
Guadeloupe. As a check of the operating conditions, 
the radar easily detected striations over Arecibo when 
the heater frequency matched the plasma frequency of 
the ionosphere, a condition that will not be satisfied 
by the SPS system (Section 4). The results are consis- 
tent with currently accepted theory (Section 4.5). 


1.3 COMMUNICATION AND NAVIGATION EFFECTS (See Section 

3.4, 5.3) 

Studies show that if striations are produced by the 
SPS system, radio frequency interference will be in- 
creased substantially in patterns that are predictable. 
Interference-free regions for ground-to-ground communi- 
cation links exist at magnetic latitudes above 60° or 
70° depending on the height of the striations. New 
England is in an interference-free region for striations 
above 250 km. Effects at lower altitudes are likely to 
be similar to sporadic-E effects except that they are 
permanent (Section 3.4). The effects produced by ther- 
mal self- focusing are limited for a 5 GW SPS system to 
angles of 5° or less between the SPS beam and the mag- 
netic field lines. This limits' the striations produced 
by thermal self- focusing to a narrow area extending from 
Miami, Florida to Brownsville, Texas. The radio fre- 
quency interference from this source will be mostly north 
south across the Gulf of Mexico.- 



1-2 


1.4 TEST PROGRSiM (See Section 5, particularly 5.9) 

An HP heater to equal or exceed the SPS radio wave 
heating (Section 5.1) can be built and operated in 
the next two years (Section 5.8) at an estimated 
cost of $4 million (Section 5.7). A crucial part 
of the 'test program is a series of diagnostic obser- 
vations to insure that the results can be scaled to 
S-band (Section 5.2, and 5.4), The test program 
should be conducted in the vicinity of the Arecibo 
Observatory (Section 5.6, 5.9). 
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2. RECOMMENDAT I ONS 


2.1 An HP heater should he designed, built and operated 
to reproduce and exceed the radio wave heating of 
the ionosphere by the 5 GW SPS system. 

2 - 2 The heating effects should be mapped and measured 
by an incoherent scatter radar, by photometers and 
by ionosondes so that the scaling from HP to S— band 
can be made with confidence. 

2.3 The effect of the ionospheric heating on communica- 
tion and navigation systems should be measured 
directly . 

2.4 The heater should be located in the vicinity of the 
2\recibo Observatory to take advantage of the exist- 
ing diagnostic capabilities and to minimize the cost 
■of the program. 

2.5 Theoretical and experimental investigations should be 
continued into the physics and phenomena of nonlinear 
wave-plasma interactions , plasma turbulence , and the 
ionospheric processes associated with substantial local 
heating. 



3-1 


3 . AHALYSIS OF SPS BEAM/IONOSPHSRE- INTERACTION 


Electromagnetic waves propagating through a plasma, are 
reflected when the wave frequency equals the local plasma 
frequency fp^/ given by 

fp^ - 9 yF . . ' (3.1) 

where n is electron density in electrons per meter cubed. 
The maximum plasma density thus defines a critical fre- 
quency, or penetration frequency, which for the ionosphere 
is usually found at the peak of the P region with typical 
values of 3 to 10 MHz, depending on time of day. Because 
microwave frequencies are much greater than ionospheric 
plasma frequencies , most of the radio wave energy passes 
through the ionosphere. The remainder goes into under- 
dense ionospheric heating. 

Ionosphere-microwave beam interactions can be divided into 
two general categories : resistive heating effects and 

self- focusing instabilities. The results of resistive 
heating by radio waves on ionospheric temperatures , elec- 
tron densities, and airglow emissions have been investi- 
gated in detail by Perkins and Roble (Appendix 1) using 
mimerical ionospheric models and heat balance codes - 
Independent theoretical derivations of self-focusing of 
radio waves in the ionosphere are given by Perkins and 
Valeo (1974), Abramovich (1977), and Cragin ^ a^. (1977). 

These theories are summarized and interpreted in terms of 
resultant physical effects for underdense ionospheric 
heating. 


3.1 RADIO WAVE HEATING OF THE IONOSPHERE 
(from Perkins and Roble, Appendix 1) 


The volume heat source Q resulting from colli sional 
(ohmic) heating of a plasma is given by 


Q 




(V 


ei 


+ V ) 
en 


(3.2) 


v;here Eq is the peak electric field amplitude, f^ is the 
radio wave operating frequency, and and are the 



electron-ion and electron-neutral momentum transfer 
collision frequencies/ respectively. The collision 
frequencies are functions of temperature and must be 
self-consistently determined ^in any reasonable theory. 

Because ohmic heating, scales as the inverse square of 
frequency, test programs designed to produce heating 
effects, comparable to a solar power satellite microwave 
beam are scaled to lower frequencies , and thus lower 
transmitted powers . Two restrictions on a scaled ex- 
periment are: (1) the ionosphere must remain under- 

dense, and (2) nonlinear resonant plasma instabilities- 
should be avoided. These parametric instabilities re- 
present a strong heat input mechanism which will not 
effect microwave transmissions. These two conditions 
are generally satisfied if we require that the radio 
wave frequency be at least twice the critical frequency 
of the ionosphere. -Because of differing plasma condi- 
tions, the results of underdense ionospheric heating 
are analysed separately for the different ionospheric 
regions . 

In the lower ionosphere, at altitudes near 100 km, -the 
electron heating is controlled by a local balance be- 
tween electron heating and electron loss processes. At 
these altitudes -thermal conduction is unimportant. A 
1 MW underdense heating experiment from an assumed 
antenna of 300 meters diameter at Arecibo produces 
heating effec-ts of 100 km comparable to the SPS micro- 
wave beam. The most dramatic heating effects occur in 
the. E region, where thermal runaway occurs. The elec- 
tron temperature will increase from 200 °K to approxi- 
mately 1000 °K and the electron density' will increase by 
a factor of about three. 

At F-region heights thermal conduction plays a dominant 
role in the heat balance equations . The key scaling 
parameter then corresponds to energy deposited along a 
magnetic flux tube per unit tim.e. Evaluated for 300 km 
altitude, this parameter equals 2.8 3nn for the SPS and 
0.84 km for Arecibo and fQ = 15 MHz. Thus it would take 
a 3 MW transmitter at 15 MHz, or a 2 MW transmitter at 
10 MHz, to simulate the effects of the SPS accurately. 
However, a 10 MHz heating experiment could operate only 
. at night, when the ionospheric plasma frequency was low. 
Although the length of heating along a magnetic field 
line for the Arecibo experiment '21 km) is much greater 
than the corresponding value of the SPS microwave heating 
5.6 km),, both cases fulfill- the essential condition 
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that this scale he small compared to the size of the 
ionosphere- F-region electron temperatures are pre- 
dicted to increase by 200-500 ®K- along magnetic field 
lines passing through the radio wave beam with enhanced 
airglow emission- also excited.. 


3.2 THERMAL RXMAWAY 

Strong radio wave heating of the ionosphere is capable 
of inducing a thermal runaway in the D and E regions 
(Holway and Meltz, 1973) . Thermal runaway can occur 
because the electron-neutral collision frequency 
increases with increasing electron temperature. As a 
result/ as the electrons get hotter, the heat input Q 
becomes larger. In addition, the rotational excitation 
of N 2 and O 2 , which Holway and Meltz used as the princi- 
pal energy loss mechanism, actually decreases with in- 
creasing electron temperature. In the Perkins and Roble 
work significant improvements have been made to the ori- 
ginal. Holway and Meltz study by including energy loss 
through vibrational excitation of N 2 and O 2 , as well as 
energy loss resulting from 0(^P) excitation. The mono- 
tonic rise of electron temperature predicted by the ori- 
ginal theory is stopped by vibrational energy losses, 
which increase rapidly for electron temperatures above 
1000 °K. These vibrational excitations then can lead to 
observable airglow emissions. 

The time scale for electron heating and cooling asso- 
ciated with a thermal runaway is a sensitive function - 
of altitude, equal to 2.5 x 10“^ sec at 100 km' and 
1 X 10“1 sec at 120 km. Threshold for achieving thermal 
runaway at 100 km is just satisfied in the center of the 
beam for a 15 MHz, 1 MW Arecibo underdense heating -ex- 
periment. Scaling to microwave frequencies, this corres- 
ponds to an effective SPS power density of approximately 
20 mW/cm2 . To assure that thermal runaway will occur in 
the experimental test program, a 2 or 3 MW Arecibo heat- 
ing experiment is essential. 


3.3 THERMAL SELF-FOCUSING 

As an electromagnetic wave propagates through a plasma 
small natural density fluctuations cause a variation in 
the index of refraction of the medi-um, resulting in a 
slight focusing and defocusing of the wave. The electric- 
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field intensity increases as the incident wave refracts 
into regions of comparatively xinderdense plasma. Ohmic 
heating then generates a themnal pressure which drives 
plasma from these focused regions , amplifying the ini- 
tial perturbation. This process is shown schematically 
in Figure 3.1. The self- focusing instability continues 
until hydrodynamic equilibrium is reached. 

-Themnal self- focusing does not have an absolure thres- 
hold in the usual sense. The focusing process, organizes 
the plasma into large-scale field-aligned irregularities 
or striationsl The threshold for creating these -stria- 
tions depends directly on the striation width Xj_. As 
the electric field intensity of the incident radio wave 
increases, the resulting striations become narrower. 

The threshold field is given by 

^ (7.2 X 10^^)T^^f ^f cos^a , (3.3) 

m e o c X 

where T^'is the electron temperature (in ®K) , f^ is the 
critical ionospheric frequency, and 9 is the angle be- 
tween the plane orthogonal to the incident wave propa- 
gation and the magnetic field. For typical values of 
the SPS microwave beam, T 3 = 1200 ®K, f^ = 2.45 GHz, and 
f( 2 - ~ 8 MHz , then 

Pth ^ (4.1 :k: 10^^)X^“^ cos^e . (3.4) 


For Arecibo, 3 = 50°,. so that 




For- P(SPS) = 
?(SPS)' = 

For 3 == 89°, 
P(SPS) = 


20 mW/cm^, 
200 mW/cm^, 

=- (1.25 
20 mW/cm^, 
200 mW/citi^, 


960 meters r 
X^ 540 meters. 

X 10 ^^) X,'^ and 
X^ ^ 158 meters; 

X 112 meters. 

X 
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The spatiaL amplification for self-focusing is exp ^) / 

where D = kj_^L/2kQ, or D 3.8 x lO'^/Xj^^' for f^ =■ 2.45 GHz 
and L 100 km, the ionospheric scale height.. Then there 
is only small spatial amplification for X_,_ > 200 meters. 
Using this approximation, for an incident power density 
of 20 mW/cm^.^. self— focusing can occur for 9 > 88.4®; 

for 200 mW/cm?-, :-we have 84,9®. This limits the inter- 
ference problem due to thermal self-focusing to SPS beam 
angles nearly parallel to the magne-tic field lines,' 

Theoretically then thermal self-focusing will be excited 
for microwave power 'densities comparable to those of the 
SPS beam and should be experimentally studied. Also, it 
is not well "known what additional plasma phenomena accom- 
pany self-focusing radio waves in .the ionosphere. 

Using this theory, what results are expected from the 
Arecibo- interim test? Defining threshold by the criteria 
for adequate spatial gain we require that D > 1. Then 
the threshold powers for' Arecibo are ; 

P(S-band: 2380 MHz) = 104 mW/cm^ 

?(430 MHz) = 10 mW/cm^ 

and ,P{HF: 8 MHz) = 1,3 x 10“^ mW/cm^, 

compared to transmitted power' densities of the limited 
test program, 

P (S~band) 

P(430 MHz) 

and P(HP) 


2 

=. 1 mW/cm 

0 

~ ,35 mW/cm" 

= 2.3x10^ mW/cra^ . 


Clearly 'then no self- focusing is ejected for 3— band 
or 430' MHz radio waves. On "the other' hand, the HF 
power density is just above threshold so that some 
effect might have 'been produced although no echoes 
were observed in Guadeloupe. 

In addition to the above conditions,, another restric- 
tion must be imposed on the incident radio wave- If 
the spatial gain condition were relaxed, the' threshold 
striation wid-th would become larger, and thus easier 
to maintain. However, the striation width cannon be- 
come .any longer than the radio wave beamwidth. At 
250 km the S-band. beamwidth is approximanaly 200 meters, 
•die 430 MHz beamwidth is approximately 600 meters, and 
the HF beamwidth is- approximately 30 km. 
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3.4 EFFECTS ON COMMDNICATION AND NAVIGATION SYSTEMS 

"Assess the effects on existing communication and navi- 
gation systems if the microwave beam exceeds the power 
density limits for nonlinear interactions with the 
ionosphere . ” 

The assessment depends on the production in the iono- 
sphere above about 110 Ton altitude of striations that 
scatter radio waves, and on irregularities inserted- 
into the earth-ionosphere waveguide below about 110 km 
that propagate low frequency waves . 'Whether the stria— 
tions are produced by nonlinear interactions or by linear 
interactions (e.g., ohmic heating) makes little differ- 
ence- Striations will scatter radio waves regardless of 
their production mechanism. Above about 110 km the 
striations have a shape ' and an orientation that limits 
their effectiveness- as scatterers and the geographic 
region in which they provide ground-to-ground links . 

The effects of the striations and of the imperfections 
in -the waveguide, are considered separately in -the 
following. 


Striations in the lower B recrion 


If produced', striations in the .E region (around 100 km 
altitude) will act like natural sporadic-E clouds and 
scatter isotropically radio waves in the frequency range 
2 to perhaps 200 MHz and perhaps higher. As wi-th spora- • 
die E, TV signals for example would be received from 
stations up to about 2000 km distant provided -the scat-- 
tarers are near the mid-path or at least visible to both 
transmitter and receiver. As with sporadic E, -the signal 
scattered decreases rapidly v/ith increasing frequency. 
Channel 2 at the low end of the VBF TV-band is most often 
seen by this' mechanism, and usually as- an interfering, 
ra-ther' 'than useful, signal. 


Striations above 110 km altitude 


If produced, striations in the ionosphere above about 
lio km will be elongated along 'the earth's magnetic 
field and -the scattering will be very aspect sensitive. 
The shape- of the striations is the result of the high 
conductivity along the field line compared to that across 
the field line. The striation may have an excess or a 
deficit of electrons compared with ambient i 
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The scattering cross-section for the striations over 
Platteville, although produced by anomalous absorption 
associated with plasma instabilities, suggests a fre- 
quency dependence that should apply to striations pro- 
duced by other means. The scattering cross-section de- 
creases relatively slowly (less than 10 db) from 40 jyiHz 
to 90 MHz> drops an additional 20 db by 157 MHz and then 
slowly (less than 10 db) to 435 MHz. An attempt to ob- 
serve the backscatter at 1300 MHz was unsuccessful. This 
behavior is shown in Figure 3.2 taken from Minkoff (1974). 

The aspect sensitivity is best described by thinking of 
specular reflection (angle of incidence equals angle of 
reflection) from a line segment where the line is along 
the magnetic field and the length of the segment deter- 
mines the angle about the specular direction in which 
significant signal, is observed. Note that the angle of 
incidence is fixed by the direction of the arriving 
wave normal and the direction of the magnetic field, 
but that the angle of reflection is satis-fied by a 
cone’ with axis on the field line and half angle equal 
to the angle of incidence. An illustration of the 
intersection of these cones with the earth’s surface 
is shown in Figure 3.3. Thus a transmitter located 
anyplace along the contour 96® (angle of incidence, 

6® from normal to the field line) will be received by 
field-aligned scatter- at 230 km over Platteville any- 
where along the contour 84® (angle of reflection, 6® 
from normal to the field line). The reception sites 
will be in a zone along contour 84° with the width of 
the zone decreasing with increasing frequency because 
the length of the line segment (the striation size 
along the field) , although fixed in physical length, 
is longer when measured in wavelengths at the higher 
frequencies. Since the scattering is occurring from 
heights around 230 km, the point-to-point path lengths 
on the ground may be as long as 2500 km.. The scatter- 
ing is described in detail by Fialer (1974) and by 
Stathacopoulos and- Barry (1974). 

The signals received are scattered not by a single 
striation but by many striations. They therefore 
introduce noise into the communication channel. The 
observed bandwidths for the Platteville experiments 
were a few kilohertz and this is consistent with the 
multipath effects associated with the volume of scat- 
terers contributing to the received signals. 'The 
narrow band (few kilohertz) character of the scatter 
process makes it of limited use for communication pur- 
poses and establishes it as a possible source of inter- 
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FIGURE 3.3 Contours of bhe earth intersection of rays 

originating at a point 250 km above Platteville I 
and meeting the geomag'netic field ah bhe angle c 
indicabed (Pialer 1974) . 
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ference for broad-band transmissions (e.g., TV, FM 
radio). The interference region, of course, is subject 
to the limits imposed by the aspect sensitivity. 

Finally, since the scattering is field-aligned, there 
will be a region of the earth free of interference 
because ' the' dip angle (the angle the field makes with 
local horisontal) of the magnetic field is large enough 
that any signals transmitted from the groiand and scat- 
tered in the .ionosphere do not return to earth but are 
scattered into space. In this case the cone described 
above on which the scattered signals may be received 
does not intersect the earth. The maximum dip angle, 

which will provide ground-to-ground communication 

links via scattering from field-aligned striations at 
an altitude h for an earth of radius r is 


. -1 
sin 


The dip angle I is related to the magnetic latitude X 
by 


tan X = -r tan I 


so 'that for scatterers at 200 or 250 km, Imax 77^1 
or 75°6 and X is 65 °3 or 62 ?8. Figure 3.4 shows the 
regions of the earth (shaded) that are free of inter- 
ference from field-aligned scatterers at 250 km alti- 
tude for ground-to-ground links. Note that it includes 
the northern U.S.A. , including all of New England. The 
lines on the figure defining the limits of reception 
are lines of constant magnetic latitude. 

A 5 GW SPS system will produce striations only when the 
angle between the SPS beam and the magnetic field is 
less than about 5® (Section 3.2). This limits stria- 
tions produced by this mechanism to low latitudes 
(roughly a line from Miami, Florida to Brownsville, Texas) 
and the radio frequency interference will be mostly be- 
tween stations north and south of the Gulf of Mexico and 
between stations ■ south of the Gulf of Mexico. 




IGURE 3.4 Region in which fxeld-alxgned scatter 
from 250 Ion altitude can reach the 


ground . 










Earth-ionosphere wavecniide 


SPS heating of the D region (50—100 km altitnde) will 
produce imperfections in the earth— ionosphere waveguide 
that will affect navigation systems (e.g.. Omega/ 

LORAN C) and any systems in the frequency bands 10 KHz 
and 100 KHz. The imperfections will essentially cast 
shadows down the guide and these shadows will introduce 
errors into the navigation systems. The errors should 
be measured both in magnitude and in the geographic 
regions affected by an imperfection. If the errors are 
stable in time and location as measured in the test 
program, a means of applying corrections should be 
available . 
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4. LIMITED TEST PROGRMI 


"Perform a limited test program at the Areciho facility 
using existing equipment- In particular the tests should 
be performed at S-band, using the maximum power and power 
density available at the facility- The test results 
should be compared with theoretical predictions as given 
in Task A. " 

4-1 SUMMARY OF LIMITED TEST OPERATION AND RESULTS 

A limited test program was conducted between 2-15 June 
1977 using the existing ionospheric heating facilities 
at the Arecibo Observatory and a portable radar placed 
to observe at normal incidence to- the earth ' s magnetic 
field lines any field-aligned irregularities of electron 
density produced over Arecibo- The Arecibo Observatory 
is operated by Cornell University with the support of the 
National Science Foundation. The portable radar was pro- 
vided and operated by the Aeronomy Laboratory of NOAA at 
Tarare near Point-a-Pitre on the French island of Guade- 
loupe made available by the Institut de Physique du Globe 
in Paris . 

The results are that the available transmitters (HF, 

430 MHz^ 2380 MHz) and the 1000-foot antenna produced no 
striations that could be detected at Guadeloupe where the 
minimum detectable cross-section was 10~3 meters squared 
(10“3m2). The radio wave heating of the ionosphere for 
the available transmitters is about 1 % for HF, 40% for 
430 MHz, and 5% for 2380 MHz of the ohmic heating that 
the proposed 5 GW SPS system will produce- The horizontal 
dimensions of the heated ionosphere are 25 km at HF, 0-7 
km at 430 MHz and 0-2 km at 2380 MHz. The latter two are 
less than the predicted scale size to produce thermal self 
focusing instabilities and the power densities are less 
than the critical power densities predicted for thermal 
self- focusing instabilities. The negative results are 
consistent with the theoretical predictions. 

As a check of the system strong echoes (cross-sections 
up to the order of 103m2) from striations over Arecibo 
' were observed at Guadeloupe when the Arecibo heating 
frequency was equal to the plasma frequency at some 
height in the ionosphere. In this case the striations 
are the result of instabilities based on resonances be- 
tween the natural plasma frequency and the heater fre- 
quency, a problem that will not be encountered with the 
SPS. 
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The results are detailed below followed by a description 
of the experimental arrangements and the interaction 
VO lime . 

4.2 RESULTS 

Heating by HF waves above the critical frequency 
Result . 

Radio wave heating by HF waves above the critical fre- 
quency of the ionosphere (the maximiim plasma frequency 
of the ionosphere) at a heating level of 1% of the ohmic 
heating predicted for a 5 GW SPS system did not produce 
field-aligned striations with radar cross-sections of 
10'“3ni2 or more. The scattering volume is approximately 
lO^^m^ -for this case so the scattering cross-section per 
unit volume is less than 10“^^m2/m^. 

When the heating frequency is more than twice the local 
plasma frequency there are no known resonance effects. 
This condition was achieved on two days, 8 and 13 June. 

On 9 June the heating frequency determined by the opera- 
ting conditions of the transmitter was well above the 
critical frequency but not above twice the critical fre- 
quency. The signals are averaged over the times given 
in Table 4.1 to improve the sensitivity of the measure- 
ment. The column headings are defined and the sources 
of the entries are given below. 

The critical frequencies are read from ionosondes located 
at the Observatory and 10 km north of the Observatory. 


Power density calculation. 

The power density at 250 km is calculated from the 
usual expressions, 

P^(r) = (4.1) ■ 

^ 4TTh^ 

where P^ is the transmitter power (60 to 100 kilowatts), 

G is the antenna gain (160 at 7 MHz, 220 at 
8 MHz, and 290 at 9 MHz), 



Date 

1977 

Time 

AST 

Heating 

Frequency 

Mtiz 

Critical 

Frequency 

MHz 

08 June 

0310-15 

9.37 

3.6 


0316-27 

9.37 

3.5 


0329-51 

9.37 

3.4 


0415-37 

9.37 

3.3 


0438-0500 

9.37 

3.1 

i 

09 June 

0320-0417 

7.80 

4. 3-4.1 


0455-0532 

6.58 

3 . 9-4 . 2 

1 3 June 

0309-37 

7.80 

3.7 


0343-0451 

6.875 

3.3 


Power 
Density 
at 250 km, 

mW/cm^ 

Ohmic Heating 
as a Fraction 
of 5 GW SPS 
Heating 

Cross-section 
Observed at 
Guadeloupe is 
less than m^ 

2.3 X 10“^ 

0.8% 

2.1 X 10"^ 

2.6 X 10“® 

1 . 0% 

2.8 X 10”^ 

2.6 X 10*“^ 

1.0% 

3.7 X 10“^ 

3.0 X 10""^ 

1.1% 

3.0 X 10"^ ’ 

3.0 X 10”^ 

1.1% 

4.0 X 10"^ 

2.2 X 10"^ 

1.0% 

1.4 X 10"^ 

1.4 X lO"*^ 

■ 0.9% 

1.0 X 10“^ 

2.9 X 10“® 

1,3% 

2.2 X 10"^ 

1.4 X 10"^ 

0-9% 

0.7 X 10“^ 
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a (taken as 1 for night observations) is the 
fraction of the power that is transmitted 
through the atmosphere, (1-Q!g^) is the radio 
wave attenuation by the atmosphere, and 

h is the altitude. 


Frequency scaling of ohmic heating. 


The scaling of the ohmic heating from the operating 
frequency of the experiment to the SPS frequency is 
derived from the expression for ohmic heating of a 
plasma, • 


Q 



(v^i + 




(4.2) 


or 


Q a P^/f^ 


(4.3) 


where E is the peak value of the electric field of the 
° heating wave, 

f and f are the operating frequency and the plasma 
frequency , and 

and are the collision frequencies of the 
electrons with the ions and with the neutrals , 

so that the ratio of ohmic heating produced by the test 
wave, Q, and the- 2450 MHz wave of the SPS is 


Q 

^SPS 


P f 

D SPS 

P 2 

DSPS f • 


(4.4) 


and the power density of the SPS 


is taken as 20 mW/cm 


2 


Guadeloupe radar sensitivity. 

The limit on the cross-section a observed at Guadeloupe 
is calculated from the radar equation and the expression 
for the minimum detectable power of the receiver. The 
ratio of these is the signal-to-noise (S/N) , the quantity 
observed at Guadeloupe. 
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S 


a P_G 
c T 


4tt 


a • 


O' A 
c 

4nr'‘ 


(4.5) 


N = kTB/n 


(4.6) 


where a 

c 


is the factor representing the transmission 
line loss taken as 0.67, 


is the transmitter power 1.44 x 10'^ watts, 

G is the antenna gain, 420 and 

A is the antenna collecting area 1200 m^ for 
an array of 4 x 40 half-wave dipoles above 
ground. 


r is the slant range to the target (660 km) 
from the radar on Guadeloupe to the iono- 
sphere directly above the Arecibo Observatory, 

k is Boltzman' s -constant 1.38 x 10“^^, 


T is the system temperature, T^, + where 

the receiver temperature Tj^ = 1000® and the 
sky temperature is 6000°, and 

B is the bandwidth of the receiver, 4000 Hz 
for a transmitted pulse length of 200 ps and 
n is a coherent integration factor that helps 
the sensitivity. Values of n used in the 
observations range from 1 to 4. 


The' observational conditions and results are summarized 
in Table 4,1. 


Heating by 430 MHz waves 

Radio wave heating of the ionosphere by 430 MHz waves 
at a heating level equal to 40% of the ohmic heating 
predicted for a 5 GW SPS system and in a volume with 
horizontal dimensions of 1 km or less rather than the 
5 km of the SPS system did not produce field-aligned 
striations with radar cross-sections of 10~3m2 or more. 
The scattering volume is approximately lol-m^ for this 
case, so the scattering cross-section per unit volume 
is less than lO^l^m^/m^. 
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The Arecibo transmitter at 430 MHz operates only in a 
pulsed mode , but will produce pulses as long as 10 ms . 
The thermal time constant at 250 km altitude is about 
2 ms. The 10-ms pulse, therefore, brings the iono- 
spheric plasma into a ’’steady-state" during the pulse. 
The observational conditions and the results are sum- 
marized in Table 4.2 using column headings defined for 
Table 4.1. 


Heating by S-band waves (2380 MHz) 

Radio wave heating of the ionosphere by 2380 MHz waves 
at a heating level equal to 5 % of the ohmic heating 
predicted for a 5 GW SPS system and in a volxmie with 
■horizontal dimensions of 0.2 km rather than the 5 km 
of the SPS system did not produce field-aligned stria- 
tions with radar cross-sections of 10~'3m2 or more. 

The Arecibo transmitter operated in a continuous wave 
mode; the power is fed into the 300-meter dish but only 
200 meters are illuminated. The resulting beam is a 
cylinder of diameter 200 meters at the heights of 
interest. The observational conditions and the results 
are siunmarized in Table 4.3, using column headings de- 
fined for Table 4.1. The scattering volume is approxi- 
mately 10l0m3 for this case, so the scattering cross- 
section per unit voltime is less than 10“12j^2/'jn3 ^ 


4 . 3 EXPERIME13TAL ARRANGEMENTS 
Arecibo 


HE facility. 

The HF transmitter built by the Office of Telecommuni- 
cation Sciences of the Department of Commerce is connected 
through six- inch coax- line to a log-periodic antenna feed 
mounted above the 300-meter dish. The transmitter pro- 
duces up to 140 KW of power at frequencies from 4 to 12 
MHz and beams it vertically in either ordinary or extra- 
ordinary modes of polarization with antenna . gains esti- 
mated at 160 for 7 MHz, 220 for 8 MHz, and 290 for 9 MHz. 
The gain increases more rapidly than frequency-squared 
because the blockage of the aperture by the feed support 
structure and cables decreases with increasing frequency. 



TABLE 4.2 


Date 

1977 

Time 

AST 

Heating 
Frequency/ 
Pulse Length 

Critical 

Frequency 

MHz 

Power 
Density 
at 250 ton, 

mW/ cm^ 

. Ohmic Heating 
as a Fraction 
of 5 GW SPS 
Heating 

Cross-section 
Observed at 
Guadeloupe is 
less than m2 

13 June 

1301-43 

430 MHz/ 

6. 5-6, 6 

0.35 

40% 

4.0 X 10“^ 



10 ms 





14 June 

1827-1908. 

430 MHz/ 

9. 5-9.0 

0.35 

40% 

1.8 X 10“^ 



10 ms 








TABLE 4 

.3 







Power 

Ohmic Heating 

Cross-section 



Heating 

Critical 

Density 

as a Fraction 

Observed at 

Date 

Time 

Frequency 

Frequency 

at 250 km. 

of 5 GW SPS 

Guadeloupe is 

1977 

AST 

MHz 

MHz 

mW/cm^ 

Heating 

less than m2 


12 June 


1044-1115 


2380 


5. 5-5. 8 


13 Juni 


1046-1115 


2380 


5. 3-5. 8 


5% 
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The HP waves are refracted by the ionosphere when the 
wave frequency is near or less than the critical fre- 
quency, but for frequencies well above the critical 
frequency the refraction effects are negligible for 
our purposes. 

The HP beamwidth is about 8° at half-power points and 
this, corresponds to horizontal dimensions of 30 kilo- 
meters for the heated volume at an altitude of 240 km. 
The vertical dimension of the heated volume is about 
10 kilometers centered near the peak of the electron 
density profile. 


430 MHz facility. 

The 430 MHz transmitter operates in a pulsed mode with 
peak powers of 2 megawatts . A waveguide connects the 
transmitter through movable joints to a line feed de- 
signed to correct for the aberration of the spherical 
reflector. The antenna gain is 10^, the beamwidth is 
0.2®, and the beam can be pointed anywhere in a 40® 
cone with vertical axis . For our purposes the beam 
was pointed almost vertically. The zenith angle was 
4° (limited by the physical presence of the HP feed) , 
the azimuth angle was magnetic north (353°). The trans- 
mitter was operated in a long (10 msec) pulse mode. The 
thermal time constant of the ionospheric plasma .is about 
2 ms. During the long pulse the plasma was heated and 
almost reached a steady-state value. The diagnostic 
observations at Arecibo and Guadeloupe were examined 
with proper allowance for heating of the plasma during 
the long pulse and the relaxation of the plasma during 
the interpulse period. 


S-band facility. 

The S-band transmitter operates in a continuous-wave 
mode with power up to 500 kW. The frequency used was 
2380 MHz. The transmitted power is delivered to a line- 
feed to correct for spherical aberration. Two hundred 
meters of the 300-meter dish are illuminated. The 
heights of interest are in the near-field of the antenna 
where the beam if a cylinder of diameter 200 meters. The 
beam was pointed almost vertically. The zenith angle 
was 5°, the azimuth angle 112.5® (the direction toward 
Guadeloupe) . 
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Guadeloupe 

Site. 

Guadeloupe was selected as the site to observe any 
striations produced over Arecibo because the geometry 
is right, i.e., at Guadeloupe the radar sees striations 
at normal' incidence to their long axes . Any striations 
in the ionosphere above about lio tan that occur, natu- 
rally or are produced artificially will be elongated along 
the direction of the earth ' s magnetic field simply be- 
cause the magnetic field constrains the charged particles 
to spiral about the field lines and it is therefore easier 
for electrons to travel along the lines than across them. 

An elongated target of charged particles will have a 
scattering coefficient that varies with the aspect angle 
wiidi which it is observed and the coefficient maximizes 
when the striation is viewed in backscatter at normal 
incidence to the long axis. This was well mapped in the 
Platteville experiments ( Radio Science 9, November 1974). 

If striations are produced in the ionosphere over Arecibo 
they will be elongated along the magnetic field line and 
observed in backscatter anywhere in the plane perpendi- 
cular to that field line. Since the field line direction 
changes slowly with height one picks a height where stria- 
tions may be expected and generates the plane perpendicular 
to the field line at that height. The map of Figure 4.1. a 
has lines representing the intersections with the earth of 
the planes perpendicular to the field lines at 110, 220 
and 330 tan over Arecibo. For striations at heights from 
220 to 250 km Guadeloupe is a good location. The site 
near Pointe-a-Pitre has a line-of- sight perpendicular to 
the earth's magnetic field over Arecibo at a height of 
230 3cm, as shown in Figure 4.1.b. 


Radar operation. 

The radar placed on Guadeloupe operated with the following 
parameters (nominal ) : 


Frequency 

Peak Power 

Pulse Lengths 

Interpulse Period 

Receiver Noise Temperature 

Estimated Sky-Noise 


50 MHz 
15 kW 
20-200 ms 
10 ras 
1000°K 
6000°K 




Arecibo 

Observatory 


Guadeloupe 

radar 


FIGURE 4.1.b 


Arecibo Guadeloupe Geometry 
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Antenna An array of 4 lines of 40 half-wave 
dipoles over the groiond. The cal- 
culated parameters are : 

Collecting area = 1200 

Antenna gain = .420 

Vertical beamwidth' - 10° (2-way) 

Horizontal beamy/idth = 2° (2-way) 

Elevation angle = 16° 

Azimuth angle = 292° 


Data. 

The data are recorded on magnetic tape and in pictures 
of the A scope (range vs, intensity) and of the Doppler 
shifted spectrum of a target. For our purposes the 
signals at the ranges of interest were integrated over 
tens of minutes to improve the sensitivity of the mea- 
surement . 

The system was checked by observing striations' produced 
over Arecibo when the heating frequency matched the iono- 
spheric plasma frequency. Field-aligned echoes over 
Arecibo with up to 30 db signal-to-noise ratios were 
observed at a slant range of 660 km from Guadeloupe. 

These echoes are from striations produced by instabili- 
ties tied to overdense heating. 


4.4 THE ■IM’ERACTION REGION 

The interaction region is defined by the intersection 
of the radar pulse volume and the modified volume. The 
modified voltime is the extension by conduction and con- 
vection of the volume in which the energy is deposited. 


Heated volume 


Energy is deposited by ohmic heating in a volume that 
is defined by the heater beamwidth horizontally and by 
the variation of the electron density vertically. The 
horizontal dimensions are 30 km, 0.7 km and 0.2 km for 
the HF, 430 MHz and S-band heaters. The vertical dimen- 
sions are about 10 km for each (see Appendix 1). The 
heated volumes are essentially vertically over Arecibo. 
The 430 MHz heated volume is displaced northward by about 
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16 km at a height of 250 km; the S-band heated volume 
is displaced toward Guadeloupe by 21 km but these dis- 
placements leave the modified voliome well within the 
radar field of view. 


Modified volume 

The energy deposited in the above volumes heats the 
electrons and they move along the magnetic field lines 
rather easily stretching out the heated volume along 
the field by a few tens of kilometers (see Appendix 1) . 
The field lines over Arecibo are at an angle of 50® to 
the horizontal, and the heated volume is distorted along 
that direction. 

The volume modified by the heating transmitters is , 
therefore, a cylinder of diameter -30, 0.7, or 0.2 km for 
HP, 430 MHz and S-band, of height about 10 km and the 
cylinder is stretched out or slid along the field line 
whose direction is 50® from horizontal in the N-S plane 
and the stretching adds about 10 km both up and down the 
field line. Figures 4.2 and 4.3 illustrate the modified 
volumes for the HP and S-band cases respectively. 


Radar field of view 


The radar observes a volume defined by the beam angles 
and the pulse length. Pulse lengths of 200 ms were used 
corresponding to a resolution length along the radar beam 
of 30 km. The horizontal beamwidth of 2® corresponds to 
a horizontal dimension of 22 km at 660 km range and the 
vertical beamwidth of 10° corresponds to a dimension nor- 
mal to beam and in the vertical plane of 110 km. The 
beamwidths represent the angles between the points on the 
antenna pattern where the power is down by a factor of 
0 . 7 from the power in the direction of the maximum (i . e . , 
these are half-power beamwidths for the antenna pattern 
squared) . The elevation angle of the beam at Guadeloupe 
is 16°, over Arecibo the elevation angle becomes 22° due 
to the curvature of the earth. The volume, therefore, 
is approximately a box centered at 230 km over Arecibo 
having dimensions 30 km along the line-of-sight (elevation 
angle 22°, azimuth 292°), and normal to the line-of-sight 
110 km in the vertical plane and 22 km horizontally. 








4-15 


Scatterincr vol-ume 

The common portion of the modified volume and the 
observed volimae is the scattering volume and has 
dimensions roughly 30 km, 0.7 km, 0.2 km {for HF, 

430 MHz-, S-band) in the directipn of the radar beam, 
about lb km vertically in the meridian plane (verti- 
cal plane- containing , the magnetic field line) and 
22 km along the_field line. The limits are set_ 
'respectively' by the heater beamwidth, by the heating 
efficiency (i.e., the electron density profile), and 
the azimuth beamwidth of the radar. The common volume 
in cubic kilometers is about 6.6 x lO^, 15.4, and 4.4 
km3 for the HP, 430 MHz, and S-band cases. 

The targets , striations , may have an aspect sensiti- 
vity that further restricts the volume since the 
elongated targets have a scattering coefficient that 
peaks at normal incidence for a backscatter radar (it 
peaks when the angle of incidence equals the angle of 
reflection for separated transmitter-receiver links ) . 


4.5 COMPARISON OF EXPERIMENT AND THEORY 

Part of the ionospheric modification experiment was 
devoted to attempts to produce field-aligned irregu.- 
larities in the underdense ionosphere using th'e 
available S-band, 430 MHz and HF facilities. The 
principal diagnostic was a 50 MHz radar stationed in 
Guadeloupe, capable of receiving direct backscatter 
from 3 meter field-aligned irregularities over Arecibo. 
Such small scale -features had been previously observed 
over Platteville for overdense heating of the' iono- 
sphere, and were believed to accompany thermal self- 
focusing striations. As predicted by theory, no signals 
were seen for underdense heating of the ionosphere by 
the subthreshold S-band and 430 MHz radio waves. Nor 
were echoes observed from HF underdense heating. 

The overdense heating experiment produced some inter- 
esting results by itself. The accepted theory for the 
generation of the small-scale field— aligned irregulari- 
ties (Perkins, 1974) required the parametric decay of 
the incident radio wave into electrostatic plasma 
oscillations. Parametric instabilities are known to 
be excited for ordinary mode (0 -mode) polarization of 
the radio wave, but not the extraordinary mode (X mode) . 
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Nevertlieless , small-scale features were observed by 
the Guadeloupe radar for both 0 and X mode HF radia- 
tion/ with charges in range of the radar signal indi-:* 
eating that the transmitted polarizations were in fact 
different since 0 and X mode waves refract differently 
in the ionosT?h«-r» .... 

The study’ ef- ’field-aligned irregularities in the over- 
dense modified ionosphere was first undertaken by 
Thome (1974) using a circular HP phased array at Platte— 
villa. His results confirmed that these striations 
formed with power densities, scale sizes, and growth 
rates close to those predicted by Perkins and Valeo 
(1974) . Development of a new diagnostic technique in 
conjunction with the recent ionospheric modification 
experiment has resulted in the first detailed observa- 
tions of individual self-focused striations, as well as 
striation maps of the entire wave-plasma interaction 
region. These results, while in comparative agreement 
with theoretical predictions , suggest that additional 
plasma effects should also be considered. 

Intense high-frequency electromagnetic radiation inci- 
dent on the overdense ionospheric plasma is Icnown to 
excite parametric instabilities, enhancing electron 
plasma oscillations observable by incoherent back- 
scatter radar. These instabilities continue to be the 
subject of intense experimental study. Of importance 
here is the fact that the strength of the enhanced ' 
plasma waves directly depends on the local intensity of 
the pump electric field.. In addition, because, of fre- 
quency and wave .number matching conditions for both the 
parametric wave-plasma interaction and the incoherent 
backscatter process,, these enhanced waves are detected 
at only one specific altitude. As a result, systematic 
measurements of the enhanced plasma line, while sweeping 
the narrow backscatter radar beam across the interaction 
region, yield a two-dimensional cross-section charac- 
teristic of the local electric field intensity. The 
electric field striation maps thus generated show self- 
focusing spatial profiles and large-scale structuring 
of the illuminated plasma. Because the direction and 
rate of the radar movement are .experimentally controlled, 
the cross-sectional dimensions of the striations are 
easily measured. Alternatively, if the radar is fixed, 
the striations follow a slow natural (E x 3) drift through 
the beam, allowing a detailed study of the small-scale 
structure within individual striations- Once the size is 
known striation velocities can be calculated from these 
drift measurements . The experimental configuration is 
shown in Figure 4.4. 
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The incoherent backscatter radar operates at 430 MHz 
with a 1/6° beamwidth, corresponding to approximately 
730 meters at a typical interaction height of 250 km. 

The HP wave beamwidth is a fiinction of the operating 
frequency, ranging from 7.2° at 7.8 MHz to 10.9° at 
5.185 MHz', or 31.5 and 47.5 km respectively at 250 km 
altitude. The HP radiation, continuous wave with ordi- 
nary polarization, is fixed at vertical incidence, 
while the 430 MHz radar can be swept in zenith angle 
at a rate of up to 1.87°/minute,- or 136 m/sec at 250 km. 
Plasma line intensities were measured over a 20 kHz 
filter bandwidth centered on the pump frequency, using 
a 500 [isec radar pulse and a pulse repetition period of 
8500 [J,seG, as determined by a 6% transmitter duty cycle. 
The data rate is not continuous; data are accumulated 
in arrays of 1024 pulses, after which there appears a 
data gap of approximately 0.4 sec in which the array 
is recorded on magnetic tape and displayed real-time. 

Results of the two basic observing strategies are shown 
in Pigures 4,5 and 4.6. The first data correspond to 
the natural drift of striations through the fixed radar 
beam. This drift gives rise to a slow regular modula- 
tion of the measured plasma line intensities, with- 
typical periods of 60 to 90 seconds. Similar fluctua- 
tions have been noted previously, without interpretation. 
The second data set, taken immediately following the 
above drift measurements, shows a series of distinct 
striations observed by sweeping the radar beam north 
across the interaction region. Typical striation dimen- 
sions are measured to be approximately 1.2 km in the 
north-south magnetic meridian plane and 1.0 km in the 
east-west perpendicular plane. Striation velocities are 
of the order of 20 m/sec to the east and smaller than 
15 m/sec in the north-south direction. North-south velo- 
city measurements are complicated by a strong spatial 
dependence of the striation width. 

These results now can be compared to overdense self- 
focusing theory. Accordingly, 

{7.6 X 10^^)T-^L“^f cos^0 . (4.7) 

th e o c j. 

Let Te=1200°K, f q = f ^ = 8 MHz, L = 100 km, 0 = 50°, 
and P = 40 pw/m^ at the center of the HP beam. Then, 


1.95 km . 


(4.8) 



Beam fixed at 3,2° north 


145736 - 150037 



I 

i-‘ 

VD 


FIGURE 4.5 Natural striations drifting through the radar beam. 
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Thus the striation width is, hy definition, about 1 km 
at' vertical incidence. The incident power density is a 
function of the angle from vertical incidence, since 
the transmitted power is not uniform over the beam, but 
can be represented by a gaussian distribution. In addi- 
tion, the north-south sweeps also vary the value of 8. 

As a result, theoretical curves allowing for this varia- 
tion of can be computed. Figure 4.7 shows experimen- 
tal data representing the width of striations (east-west 
and north-south components shown separately) versus their 
angle from vertical incidence. The solid curves are 
theoretical expectations accepting reasonable values for 
the striation width at the center of the beam. Despite 
the large scatter in the data, the theoretical curves 
approximate the data behavior fairly well. 

An unexpected feature of the observations is the strong 
tendency for striations to appear in compact bunches , 
rather than occur imiformly distributed over the heated 
region. The fluctuations of signal strength during 
natural drifts of the striations through the radar beam 
suggest the presence of small scale irregularities such 
as those seen by the Guadeloupe radar. However, the 
diagnostic radar beam is too wide to directly resolve 
discrete striations of this size. 

Although the present technique was developed for the 
observation of overdense self- focusing effects, a 
similar experimental program measuring the altitude 
changes of naturally enhanced plasma waves should work 
to make comparable observations of underdense self- 
focusing striations possible. 



STRIATION WIDTH (km) 



FIGURE 4.7 Variation of striatioh width as a function 
of power and 0 - experiment and theory. 
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5. THE TEST PROGRAM 

"Define a test program for detailed, ground testing to 
fully determine the exact density levels at which the 
interactions occur, the extent of possible disruptions 
to navigation and communication systems, and any 
environmental impacts. This test program definition’ 
should include (1) measurement parameters, (2) evalua- 
tion of existing facilities at Arecibo and Platteville 
and the modifications needed at each of these sites to 
accommodate S-band testing at power densities up to 
50 to 100 milliwatts/cm2 , (4) costing information for 

such a test program to fully determine the microwave 
beam ionospheric interactions, and (5) recommendations 
on the options given. " 

A test program is defined to simulate the interaction 
of the SPS beam with the ionosphere to measure the 
effects of the interaction and to interpret the results. 
The program includes (1) an HF heater capable of genera- 
ting and delivering 15 MHz waves that produce more ohmic 
heating at all levels of the ionosphere than a 5 GW SPS 
and with a beam diameter that equals or exceeds the SPS 
beam diameter at all heights where it is critical, (2) 
diagnostics that map the changes in electron density 
and temperature and measure the intensity of any plasma 
instabilities produced, (3) field measurements to assess 
the effects on communications and navigation, (4) analy- 
sis and interpretation of the measurements, (5) associated 
studies, (6) a comparison of Arecibo and Platteville 
facilities, (7) budgetary estimates, (8) schedule, and 
(9) reco3mnended test. Each of the nine elements of the 
test program is described below. 


5.1 HF HEATER 

A simulation of the SPS beam interaction with the iono- 
sphere is necessary to keep the experiment manageable. 

The frequency scaling law for ohmic heating is known 
(Equation 4.4) and this provides a way of introducing 
the same quantity of heat into the medium from a source 
of a few megawatts rather than a few gigawatts. The 
scaling of the radio wave refraction in an ionized medium 
is also known and although this cannot, be scaled in the 
experiment the proposed diagnostics permit the scaling to 
be made in the analysis and interpretation. 
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Heater frecruency 

The frequency scaling law for ohmic heating {Equation 
4.4) shows that for the same heating (Q = QsPS^ 
test power density P that must he delivered to the 
ionosphere is 

^2 

“ ^DSPS T2 (5.1 

^SPS 

and this P^ is minimized for a given SPS system 
(PdSPS and fspg) by making f as small as possible. 

A lower limit on the frequency is set by noting that 
one does not wish to excite any parametric instabili- 
ties in the ionospheric plasma in the test, for these 
are known to produce field-aligned striations. The 
SPS system frequency is far removed from the condition 
for parametric instabilities and the test frequency 
must, therefore, also be removed from that condition. 
To satisfy this condition one requires that the test 
frequency be at least twice the highest local plasma 
frequency, the critical frequency- Typical critical 
frequencies are 3 to 10 MHz depending on time of day 
and other factors. 

We select 15 MHz as the nominal frequency that will 
satisfy the above conditions and permit operations 
for more than 50% of the time. 


Power density 

The 5 GW SPS system has typical power densities in the 
beam of 20-30 mW/cm^. The contract calls for a test 
program with power densities equivalent to 50-100 mW/cm^ . 
Meltz (Space-Based Solar Power Conversion and Delivery 
Systems Study, Vol. Ill, Microwave Power Transmission 
Studies, Raytheon Co., March 1977) expects that non- 
linear effects may occur at 20mwAm2 and above and rec- 
commends a test at these levels. Our studies (Section 
3-2) show that theimial aninaway in the E- region is likely 
to occur at 20 mW/cm2 . Perkins and Roble (Appendix 1) 
show that 20 mW/cm2 will produce fractional changes of 
electron temperature of 30% or more in the ionospheric 
F- region. Although we have no firm prediction that this 
will result in an instability, the heating clearly needs 
to be produced and watched for instabilities. 
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Considering all of the inputs given above we use as 
criteria a test power density in excess of 25 mW/cm^ 
at 300 km, 50 mw/cm^ or more at 220 km and below, and 
in excess of 100 mw/cm2 at 100 km. 


Transmitter and antenna 


The SPS beam contains about 98% of the power in a 
circle of 10 km diameter, and about 63% of the power 
in a circle of 5 km diameter. Since instabilities 
have a power threshold that must be exceeded to set 
them off, we concentrate on the 5 km diameter beam 
size where the SPS power densities are highest. The 
power density criteria given above can be produced in 
a circle of 5 km (or more) by transmitters with a 
power of 2.4 MW and an antenna of 12 modules, each 
consisting of an array of 16 x 16 crossed half-wave 
dipoles above a ground screen. The modules are 
arranged in the form of a cross as shown. Each 
module is fed 0.2 MW power. 



1 

2 


3 

4 

5 

9 

6 

7 

8 

10 


12 

11 



Antenna Modules 


For D and E region tests, mtodules 1 through 8 are 
excited producing the equivalent of 138 mW/cm^ at 
100 km in a beam just under 5 km diameter. The 
beam diameter is 5 km at 120 km increasing linearly 
with height. The equivalent power density is 50 mW/ 
cm^ at 166 km decreasing with height squared. 

For F region studies all 12 modules are excited 
producing the equivalent of 61 mW/cm2 at 200 km, 

50 mW/cm^ at 220 km, and 27 mW/cm^ at 300 km in a 
beam size that is 5 km diameter at 160 km increasing 
linearly' with height. The heater characteristics are 
simmarized in Table 5.1. 
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TABLE 5.1 
HP Heater 

(Frequency 15 MHz; Wavelength 20 m; 
Transmitter Power 2.4 MW) 



D and E Regions 

P Region 

Power 

1.8 MW 

2.4 MW 

Niaiiber of Antenna 



Modules 

8 

12 

Approximate Antenna 

24X X 24X 

32X X 32X 

Size (see Figure) 

480m X 480m 

640m X 640m 

Equivalent SPS . 



Power Density and 



Beam Diameter at 

2 


Altitude of 80 km 

216 mW/cm /3.3 km 


100 km 

138 mW/cm^/4.2 km 


120 km 

96 mW/cm^/5.0 km 


160 km 

50 mW/cm^/6.7 km 

96 mW/cm^/5.0 km 

200 km 


61 mW/cm^/6.2 km 

220 km 


50 mW/cm^/6.9 km 

300 km 


27 raW/cm^/9.4 km 
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5.2 DIAGNOSTICS 

To observe the interaction of the simulated SPS microwave 
beam with the ionosphere requires an incoherent scatter 
(IS) radar capable of mapping throughout the heated volume 
changes’ in electron temperature and electron density with 
spatial resolution of 0.5 km^ temperature resolution of 
5 %, electron density resolution of 1 % and time resolution 
better than a minute. The IS radar must be capable of 
detecting any plasma lines excited by the heater down to 
the level of the receiver noise for a receiver front end 
temperature of 100 °K. Photometers to search for changes 
in the excitation of optical line radiation and iono- 
sondes to follow the gross changes in the ionosphere are 
required. 


5.3 FIELD MEASUREMEOTS OF EFFECTS ON COMMUNICATION 
AND NAVIGATION SYSTEMS 

To establish the effects, if any, on communica'tion and 
navigation systems, we propose. measurements at field 
sites chosen to emphasize the effects (see Section 3) 
associated with field-aligned striations and earth- 
ionosphere waveguide irregularities . 

A 50 MHz radar of the type used successfully in the 
limited test on Guadeloupe (Section 4.3) is appropriate 
to look for field-aligned striations . Any radar that 
equals or exceeds in performance that used on Guadeloupe 
is acceptable and any site may be used that meets the 
aspect sensitivity requirements. Searches from two 
sites, however, are required. One to cover striations 
in the F region, the other to look for echoes in the 
E region. The E-region search should be conducted from 
a site where the aspect sensitivity from striations at 
120 km is satisfied. A radar at this location should 
look for field-aligned echoes from 120 km and non-field- 
aligned echoes from lower heights . 

To check for navigation problems , Loran-C and Omega 
receivers should be set up at fixed locations along the 
great circle from the navigation transnutter through the 
heater to look for changes in the system indicators with 
and without the heater operating. 
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The field measurements outlined here are first checks . 

If the results are negative then a wide range of possi- 
ble problems are eliminated. If any of the results are 
positive the diagnostic observations must be analysed 
to assure a proper scaling to the SPS system. Additional 
field tests may be needed to define the scope of the 
difficulties. . , ' 


5-4 AliJALYSIS AMD INTERPRETATION ‘ 

The diagnostic measurements and the field measurements 
will be analysed to define the physical effects induced . 
in the ionospheric plasma by the heater and relate them 
to the field measurements, either positive or negative. 

Mapping of the electron temperature and density with and 
without the heater will yield plasma motions and scales , 
the intensities of any irregularities, the rates at which 
changes occur, and the time constants ' for achieving a new 
"steady- state. " 

While the scaling of the heating follows that of the SPS 
beam, the refraction is not scaled and the mapping will 
be crucial in establishing that the effects achieved in 
the simulation will also occur in the SPS beam. 

Photometer observations will be useful indicators of 
both ’the usual thermal effects and the predicted thermal 
runaway in the D and E regions. 

The interpretation will be based on the results of the 
field measurements in search of communication and navi- 
gation problems and on the analysis of the diagnostic 
data to define the areas that can be eliminated from 
further study, those that may have serious impact, and 
those that require additional testing. 


5.5 ASSOCIATED. STUDIES 

Although theoretical plasma physics has been very suc- 
cessful at explaining ionospheric modification experi- 
mental measurements, it has been comparatively poor at 
predicting new results. The elegant mathematical tech- 
niques of linear plasma theoary are -relatively ‘unprepared 
to explore the jungle of nonlinear phenomena found in 
experimental observations. As a result, there is a large 
inherent uncertainty- in any discussion of future research 
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needs. The following section identifies several 
potentially important plasma effects which deserve 
further consideration. 


Ponderomotive force 

The ponderomotive forcfe is a nonlinear plasma force 
deriving from the second order particle equation of 
motion, directly dependent on the gradient of the 
electric field strength. Although the ponderomotive 
force is generally much weaker than the ohmic heating 
thermal forces associated with the microwave beam, 
sharp gradients in electric field strength accompany- 
ing self- focused striations may introduce additional 
nonlinear effects . This gradient would be emphasized 
near the striation boundaries, suggesting that pondero- 
motive forces may produce irregularities on smaller 
scales than thermal self-focusing. 

In the past instabilities driven by the ponderomotive 
force have been the most difficult plasma effects to 
theoretically anticipate. Because of the high power 
densities involved, nonlinear microwave-ionosphere, 
interactions -can generate completely unexpected plasma 
phenomena. For this reason, continued experimental and 
theoretical research of nonlinear plasma behavior is a 
necessity. To improve our imderstanding of self- focus- 
ing waves in plasmas , a self-consistent theory combin- 
ing both thermal and ponderomotive forces needs to be 
developed. 

Stimulated scatter 


Self- focusing of the microwave beam, as discussed in 
Section 3.3, is a foirward scattering mode of stimulated 
scatter driven by thermal forces. As such, only a small 
fraction of a percent of the incident radiation is 
actually absorbed or scattered away by the plasma. A 
much more serious problem would occur if the radiation, 
suffered appreciable sidescatter or backscatter. In 
this case, a significant loss of transmission efficiency 
would result. 

Stimulated scatter of electromagnetic radiation has been 
investigated as it applies to scatter from ion acoustic 
waves (stimulated Brillouin scatter) and electron plasma 
waves (stimulated Raman scatter). Fortunately, in both 
cases stimulated sidescatter and backscatter can be 
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neglected because of their much higher thresholds . 
Present theoiy states that these instabilities will 
be unimportant for the microwave power beam. However, 
as has been mentioned, theory is much better at ex- 
plaining results than predicting observations. The 
possibility of additional stimulated scatter should 
not be forgotten, in case new studies significantly 
change our ^lnderstanding of the scattering process . 


Additional effects 

The effect of microwave radiation on atypical iono- 
spheric conditions needs to be studied. The amplifi- 
cation of existing ionospheric irregularities or an 
increased tendency for unusual ionospheric conditions 
represents a serious potential problem in ionosphere- 
microwave beam interactions. In addition, the possi- 
bility for generating or amplifying travelling iono- 
spheric disturbances needs to be considered. All of 
the above effects can result in ionospheric irregu- 
larities over areas much larger than that intersected 
by the' microwave beam. 

Also, some basic physics questions remain unanswered. 

We have not discussed the possibility of radio fre- 
quency interference due to harmonic generation. The 
basic chemistry of the ionosphere could change dras- 
tically due to the large local input of thermal energy. 
Results from the most recent ionospheric modification 
experiment have left us without a viable explanation 
for the generation of small-scale field-aligned irre- 
gularities. We need to iinderstand how plasma motions 
are affected by turbulent decay processes. 

Parametric instabilities in the ionosphere have been 
suppressed successfully via frequency, phase, and 
amplitude modulation of the incident pump radiation. 
Although the spatial instabilities considered in this 
report differ considerably from parametric instabilities, 
continued research into the details of wave-plasma in- 
teractions offers interesting opportunities . 


Laboratory plasma experiments 

Detailed studies of intense ohmic heating effects in 
a plasma have been proposed based on ionospheric 'modi- 
fication research. In all cases, these proposals have 
been developed around a frequency scaled experiment. 
That is, by lowering the operating frequency an equi- 
valent amount of ohmic heating is achieved for much 
lower transmitter powers compared to the SPS microwave 
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system. Because the ionosphere is an excellent 
labor atory-without-wal Is for investigating large-scale 
plasma irregularities , the frequency scaled experiments 
offer promising research opportunities. However, addi- 
tional experimental configurations were also considered. 

A study was made of the feasibility for scaled labora- 
tory studies of ionosphere-microwave interactions using 
Johnson Space Center ' s large thermal- vacuum test facility 
(Chamber A) . An experiment was considered using micro- 
wave frequencies while scaling the radiated power and 
plasma density. Consequently a study was made to iden- 
tify the principal problem areas for such an experimental 
program, the results of which appear in a report given in 
Appendix 2 . The questions which arose out of this report 
were then studied, with the conclusion that a scaled 
laboratory experiment was not practical at this time. 

Laboratory experiments at microwave frequencies lack 
adequate spatial amplification for self- focusing to 
occur. The theoretical condition for forward scattering 
requires > Xq/s » e small. In addition, the condition 
for spatial amplification can be expressed as Xj_ < N Xq, 
where N is the number of wavelengths Xj, which go into the 
plasma scale length L, or X, = L/N. These restrictions 
then yield a relationship defining a key scaling factor, 

X^L s . (5.2) 

Thus for microwave frequencies = 2.2 GHz, Xq = 13.6 cm) 

and small angle scattering (e = 10-1 is still fairly large 
angle scattering) the plasma scale length for self-focusing 
spatial amplification is L s 13.6 meters. This large 
dimension suggests that scaled laboratory experiments would 
be only minimally useful since the working plasma scale 
length of JSC ' s Chamber A is approximately 10 meters . 

It is interesting to note that by increasing the operating 
frequency, the plasma scale length for self- focusing can 
be reduced. However, a test program using frequencies 
higher than those proposed for the -SPS microwave beam is 
not necessary, since the equivalent phenomenon of beam 
filamentation has already been studied using high power 
densities and laser radiation. These results are not 
directly applicable to the study of ionosphere-microwave 
beam interactions . 
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Although laboratory experimentation does not appear 
promising for the investigation of thermal self- focusing 
instabilities, it continues to play an important role in 
the overall study of the ionosphere. Laboratory work 
provides valuable information on the ionospheric chemistry, 
collisiohal processes, and plasma instabilities, all of 
which are important if we hope to understand how the iono- 
sphere behaves when strongly heated. Attention should be 
quickly given to any new laboratory research results which 
affect our present understanding of ionospheric processes . 


5.6 COMPARISON OF ARECIBO AND PLATTEVILLE FACILITIES' 

The requirements for the test program can be listed under 
three main headings: the heater, the diagnostics, the 

field measurements. The general requirements for each 
item and the present availability at each site is indicated 
in the table below. 


TABLE 5.2 


Heater 

Transmitter (15 MHz) 
15 MHz Antenna 
Transmission Line 
Site 

Diagnostics 

Incoherent scatter 
radar 

Photometer 

lonosondes 


Availability at 

Requirement Arecibo Platte vi lie 


2.4 MW 

640m X 640m 

~ 6 km 

120 acre flat 
rentable 


140 Kw 
None 
None 

Available 


2.0 MW 
None 
3 km 

Available 


Equivalent AO Available 


None 


Available Available 
Available Available 


Field Measurements 
Sites 
Equipment 


Available Available 
Available Available 
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The advantage is with Platteville with respect to the 
heater and the value of the advantage is about $2,5 
million. The advantage is with Arecibo with respect 
to the diagnostics and the dollar value of the advantage 
is well in excess of $10 million. Without attempting 
to price-out an IS radar at Platteville equivalent to 
Arecibo, one notes that the cost of the Arecibo facility 
was $9 million in 1963. 

There is no apparent advantage with respect to field 
measurements from the table but it should be noted that 
for a Platteville heater all of the interesting area is 
land and in the U.S.A, For an Arecibo heater much of 
the interesting area is sea with islands, few of which 
are U.S. territory. 

The advantage with respect to cost is heavily in favor 
of Arecibo. 

5 . 7 ESTIMATED BUDGET 

The estimates tabulated below are based on inputs from 
Continental Electronics , the Aeronomy Laboratory of 
NOAA, L. M. LaLonde, and a proposal from Cornell University 
to the National Science Foundation for an HF Ionospheric 
Modification Facility dated November 1975 written by W. E. 
Gordon and L. M. LaLonde with the aid of AO staff members. 

Heater in Place 

Transmitter 15 MHz, 2.4 MW 

Antenna 12 modules, 

16 X 16 crossed dipoles 

Transmission line and dividers 

Primary power 4000 kW 

Site development and 
transmitter shelter 


Field Sites and Equipment 

Operations , Associated Studies 
and Analysis 

Staff 

Site rentals , field expenses 
and travel 

Power 

Materials and supplies 

1,000 K 


$1,000 K 

750 K 
300 K 
600 K 

300 K 

$2,950 K 

$ 50 K 50 K 

$ 600 K 

90 K 
130 K 

180 K 


TOTAL ESTIMATED COST 


$4,000 X 
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BtiDGET NOTES: 

(1) The transmitter estimate depends on securing 12 
FRT-86 Navy units for modification. The units 
are presently available. 

(2) The antenna modules are 16 lines x 16 lines of 
16 dipoles of the kind described by Balsley and 
Ecklund (1972). 

(3) The transmission lines are fabricated locally 
from standard irrigation pipes. 

(4) The equipment at the field sites will be on loan. 


5.8 SCHEDULE 

The schedule outlined in Table 5.3 calls for final design 
of the heater by early 1978, the procurement and installa- 
tion of the heater by the end of 1978, tests in the first 
half of 1979 with the second half of 1979 reserved for 
further tests if warranted. 


5 . 9 RECOMMENDATIONS 

(1) An HF heater should be designed, built and operated 
to reproduce and exceed the radio wave heating of 
the ionosphere of the 5 GW SPS system. 

(2) The heating effects should be mapped and measured 
by an incoherent scatter radar, by photometers and 
by ionosondes so that the scaling from HF to S-band 
can be made with confidence . 

(3) The effect of the ionospheric heating on communica- 
tion and navigation systems should be measured 
directly. 

(4) The heater should be located in the vicinity of the 
Arecibo Observatory to take advantage of the exist- 
ing diagnostic capabilities and to minimize the cost 
of the program. 
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TABLE 5.3 
SCHEDTJLE 
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APPEHDIX 1 

IONOSPHERIC HEATING BY RADIOWAVES: 

PREDICTIONS FOR ARECIBO AND THE SATELLITE POWER STATION 

F. W. Perkins* and R. G. Roble ^ 

National Center for Atmospheric Research 
?. 0. Box 3000 
Boulder, Colorado 80307 

ABSTRACT 

The effect of resistive heating by radiowaves on ionospheric 
temperatures, electron densities, and airglow emissions is examined by 
using nximerical ionospheric structure and heat balance codes. Two cases 
are studied: (1) a 3 GHz, 10 GW microwave beam from a proposed satel- 

lite power station and (2) 1 MW and 3 MW beams of 15 MHz radiowaves 
launched by the Arecibo antenna. -By intent, these two cases have simi- 
lar intensities and geometries of resistive heating. The most dramatic 
heating effects are predicted to occur in the E-region, where a thermal 
runaway will take place. The E-region electron temperature will increase 
from 200°K to roughly 1000°K, and the E-region electron density will in- 
crease by a factor of about three. In the F-region, where thermal con- 
ductivity plays an important role, temperature increases of 200-500'‘K 
will appear along magnetic field lines passing through the radiowave 
beams . Enhanced emissions in airglow and molecular infrared lines will 
also occur. Radiowave heating, when combined with the diagnostic cap- 
abilities of the Arecibo incoherent scatter radar, will generate new 
opportunities to measure the rates of atomic physics processes and neutral 
atmosphere temperatures and composition at D- and E-region altitudes, 

*0n leave from Plasma Physics Laboratory, P.O. Box 451, Princeton, 

N.J. 08540. 

^The National Center for Atmospheric Research is sponsored by -the 
National Science Foundation. 
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INTRODUCTION 


Underdense ionospheric heating occurs when the earth's ionosphere 
is heated by a beam of high power radiowaves and most of the radiowave 
energy propagates through the ionosphere because the ionospheric plasma 
frequency remains well below the radiowave frequency. Two systems which 
have this property are considered in detail: (1) a satellite power 

station [Brown, 1973; Glasser, 1977] which transmits from geosynchronous 
orbit to ground a radiowave beam with a diameter of roughly 8 km and 
which has a power- of 10 MW at a frequency of 3 GHz, and (2) an iono- 
spheric" heating facility using the Arecibo antenna and having a power' 
of 1-3 MW with a frequency in the range of 5-30 MHz. Although the 
F-region heating of the ionosphere is generally overdense at the lower 
end of these frequency ranges, heating in the E-region is still under- 
dense and, as is shown later, it is in the E-region that the most 
dramatic heating effects take place. 

We examine the physics that occurs when a radidwave beam is incident 
on an underdense ionosphere. There are two general classes of resulting 
physical effects: resistive heating and self-focusing instabilities. 

This paper concentrates "solely on the resistive heating effects. A 
subsequent paper will develop existing self-focusing theories [Perkins 
and Valeo, 1974] in the context of underdense ionospheric heating 
experiments. 

We show that, in general, the heating effects of a satellite power 
station can be simulated by an appropriately scaled Arecibo heating 
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facility. In. addition to enabling simulations that have engineering 
interest with regard to the satellite power station, underdense heating 
by the Arecibo heating facility provides a tool for performing active 
experiments in ionospheric heat balance, optical and infrared airglow 
kinetics, and diagnostics of the neutral atmosphere above about 80 km. 
These exciting possibilities derive from two facts: the radiowave- 

induced heating is as strong as the natural heating from the sun, and 
it can be temporally controlled, allowing measurements of time constants, 
coherent detection techniques,, etc. 
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RADIOWAVE HEATING OF THE IONOSPHERIC PLASMA 


Eadiowaves propagating through the ionosphere are attenuated by 
electron-ion and electron-neutral collisions. The voliome heat source 
Q which results from these processes is 


Q = 


2 2 

0 pe 

8tt 2 
V 


(V . + V ) 
ei en 


( 1 ) 


where Eq is the peak electric amplitude, is the ionospheric plasma 


pe- 


frequency, v is the radiowave frequency, and V , and V are the electron- 

ex en 

ion and electron-neutral momentum transfer collision frequencies, res- 
pectively. These collision frequencies are functions of temperature and 
they must be self-consistently determined in any code that is used to 
calculate ionospheric effects. The electric field Eq is determined by 
the extent and power of the radiowave beam and must be evaluated for 
beams from both a satellite power station and a Arecibo heating facility. 

Satellite Power Station . We assume that a satellite power station 
will beam to- earth a radiowave that has a Gaussian shape such that 95% 
of the power is contained in a circle 8,6 km in diameter. Figure 1 por- 
trays the geometry of radiowave propagation from geos 3 rnchronous orbit 
to earth. The power is assumed to be 10 GW and the radiowave frequency 
to be 3 GHz [Erown, 1973; Glasser, 1977], The shape of the radiowave 
beam has a power flux given by 


F - c^ expC-r^/r^) 


( 2 ) 



-2 

whare =» 500 Wa and = 2.5 km. The volume ionospheric heating 
rate associated with this flux is 


Fv 


Q = 


pe 


cv 






fS_n (V .. -r V ) 
0 e ex en 


where 


AiTe^F^ 

T" 

nuj c 


2.95x10 ergs, 


Fq * 10^ ergs cm"'^ « 1000 W/m^, 

■f = F/Fq = 0.5 expC-r^/r^), and 


( 3 ) 


n 


e 


electron number density. 


The altitude dependence of Q is deterinined by how the spatial varia- 
tion of the beam is mapped into a height variation along the magnetic 
field line. We consider two arbitrary ground stations, in Boulder, 
Colorado, and Houston, Texas, whose collectors receive radiowave energy 
beamed from a satellite in geosj^nchronous orbit. The geometry giving 
the intersection of the radiowave beam with each ground station and geo- 
magnetic field line is illustrated in Figure 1, The relevant parameters 
for Boulder are D = 70“’ and S= 26.4°; for Houston they are D = 49° and 
0=9°. 


From the geometry of Figure 1, the parameter f becomes 
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where 5, = sin D/ sin 3 » z is altitude and is the altitude at 
which the center of the radiowave beam intersects the geomagnetic field 
line. The value of S, is 5,6 km for Boulder and 15.9 km for Houston. 

Arecibo Heating Facility . The resistive heating generated by an 
Arecibo heating facility can be calciilated by assuming a beam pattern 
such as that illustrated in Figure 2 and represented by the formula 


F = (-0^/0 q) 

7T0OZ 


(5) 


where P is the heating power and z is the altitude. The angle 

0Q-= ’2>/T/ka, where a is the diameter of' the antenna and k'is the wave- 

number. Formula (5) assumes the effective area of the Arecibo antenna 

2 4 2 

to be A^^^ =s 0.5 A, with A being the geometrical area, A = TTa « 7x10 m 
and a = 150 m. The volume heating rate due to radiowave heating by the 
Arecibo facility is 


Q = fe-n (V + V .) 

0 e en ex 


( 6 ) 


—16 

where- Sq = 2.95x10 ergs and 


f = 0.09 ( ) exp [-(z-z ) 2/^2] 

\ 7x10 MW m / \ ^ 


where PA denotes the power aperture product appropriate to the antenna 
and transmitter and is the altitude at which the radiowave beam inter- 
sects the geomagnetic field line. The length I is given by 
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I = 


tan D 


7 km 


( \ / 30 MH2 \ 

\ 200 km / \ V / 


(7) 


where v is the radiowave frequency and the value D = 50” has been xised 
for Areciho, 

Scaling Analysis . From the heating relations derived above it can 
be seen that a 1 MW underdense ionospheric heating experiment at Arecibo 
can produce heating effects comparable to those from a satellite power 
station, especially at altitudes on the order of 100 km. In the E~ and 
D-regions of the ionosphere, the electron heating is- controlled by a 
local balance between electron heating and electron loss processes. 

Thermal conduction can be neglected at these altitudes. At 100 km, the 
Arecibo simiilation gives f = 0.36, compared to f = 0.50 for the satellite 
power station, and with the frequency of the modifier chosen to be 
V = 15 MHz, the spatial scale is S. = 7 km for Arecibo and 5. = 5.6 km for 
the satellite power station beaming radiowave energy to Boulder. 

At F~region heights, where thermal conduction plays the dominant 
role, the key scaling parameter is fi, which physically corresponds to 
the energy deposited on a magnetic flux tube per unit time. Using the 
expressions given above evaluated at 300 km, the satellite power station 
yields f£ = 2.8 km, and for the V = 15 MHz Arecibo simulation fA =* 0.84 km. 
Thus, the 1 MW Arecibo heating is roughly one-third as intense as that 
from the satellite power station. Hence it would take a 3 MW heater, or 
a _2 MH heater at V - 10 MHz, to simulate the effects of the satellite 
power station accurately. An underdense 10 MHz heating experiment will 
work only at night, when the ionospheric plasma frequencies are small 
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enough that refraction is not a problem. Even, though the length of 
heating along a magnetic line for the Arecibo experiment ~ 21 km) 
is much greater than the value for the satellite power station 
(i =5 5.6 km) , both cases fulfill the essential condition - that this 
length be small compared to the size of the ionosphere. Finally, if 

4 

one were able to use the full Arecibo aperture = A = 7.10 m), 

the E-region heat input would double and F-region heating would increase 
by a factor of 'Pi [due to the factor 9^ in (7)]. 
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THEBHAL RUNAWAY IN THE D- AND E-REGIONS 


Holway and Meltz [1973] first showed that strong radiowave heating 
could induce a thermal runaway in the D- and E-regions. The reason for 
this thermal, runaway - is that the electron— neutron collision frequency 

increases with increasing electron temperature; Che hotter the elec- 
trons get, the larger the volume heat source Q becomes [see (1)].- More- 
over, the rotational excitation of N^ and 0^, which Holway and Meltz 
assumed to be the principal electron energy loss mechanism, actually 
decreases with increasing electron temperature thus contributing to the 
thermal" rmaway. We’ shall show in- the’ following ’’discussion’ that ‘if the 
resistive radiowave heating is sufficiently large, this model has no 
steady-state solution to the electron heat balance equation and the 
electron temperature will monotonically increase in time. 

But the- Holway and Meltz model has only limited physics ;~it speci- 
fically neglects electron energy loss via vibrational excitation of N^ 

3 

and 0^, as well as energy loss resulting from 0( P) excitation. The 
vibrational excitation processes provide an energy loss for electrons 
that increases rapidly as the electron temperature rises above 1000‘’X. 
Hence vibrational energy loss will stop the monotonic rise in electron 
temperature. In this nonlinear stage, the energy input through resistive 
heating will appear as vibrational excitation of N^ snd O^. The litera- 
ture then suggests that the energy in the vibrational excitations will 
be resonantly exchanged to vibrational modes of CO^ and H^O, and radiated 
as 4.3 ]i and 6.5 U photons [Houghton, 1967; Kumer and James,. 1974; Kumer, 
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1977a, 1977b]. Energy loss resulting from 0( P) excitations proceeds 
at all relevant electron temperatures and combines with rotational 
losses to determine the threshold for initiating the electron thermal 
runaway. 

Let us now develop the equations which govern the electron- thermal 
runaway. If we neglect thermal conduction, the electron energy equation 
becomes 
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where is the electron temperature in, energy units, is the neutral 

gas temperature, t is time, n^ is the electron density, e the electronic 

charge, m the mass of an electron, c the speed of light, to the radiowave 

frequency, F the power flux, V the electron-neutral collision fre- 

en 

quency, v'' ^ the electron-rotational excitation energy loss frequency, 
and the oxygen fine structure energy loss frequency. According to 

data, when n(0) « n(N 2 ) and nCO^), we can fit these quantities suffi- 
ciently well by the formulas [Massey ^ al. , 1969; Hoegy, 1976], 
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These formulas give the collision frequencies that are appropriate when 
= T^, The explicit dependence of collision frequencies on electron 
temperature appears in (8) , which one can cast into the nondimensioual 
form 
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Thermal runaway occurs when steady-state solution to (12) is not possi- 
ble because of the high value of A. This condition can be evaluated 
by solving the two equations 
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where C = y . These conditions find the value of A for which the 

c 

curve A^5^/(l+aO is tangent to 5^-1. The solution is 
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shown, in Figure 3. 
asymptotic formulas for are 
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For standard abundances of the molecular species and n(0) « n(N^) + n(02) 
one can transform (13) and (14) into 
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where 


"0 n(N.2) + n(Q2) 

We have used (5) and (6), appropriate for the Arecibo experiment, to 
evaluate F appearing in (13) . 

The U.S. Standard Atmosphere , [NASA et al. , 1976] has Oq = 0.043 

at 100 km yielding a = .21 for T = 200°K, Hence these numbers suggest 

n 

that a 1 MW Arecibo underdense heating experiment would be just on the 
margin of achieving a thermal runaway at 100 km and then only in the 
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center of the beam. At lower altitudes, where both Zq and the atomic 
oxygen concentration are smaller, thermal riinaway will occur. Above 
100 km, the cooling of atomic oxygen prevents a large thermal runaway. 
Clearly, to assure that thermal runaway will occur, a 2 or 3 MW Arecibo 
heating experiment is essential. 

The time scale for heating and cooling is essentially 



(23) 


The rapid heating and cooling at 100 km suggest that a pulsed experi- 
ment might well suffice to get high peak powers and thus not be as 
demanding on the average power. For example, a 10 ms pulse on the 
430 MHz. Arecibo radar could, achieve steady-state heating conditions 
and still satisfy the thermal conductivity requirements of (26). The 
problem of how to use the 430 MHz radar both as a heater and as a diag- 
nostic tool can be solved by employing separate but neighboring fre- 
quencies . 

In the collision-dominated regime, the Arecibo Thomson scatter 

2 2 

spectrum is Lorentzian, P(w) = (const)/ (to +y ) where Y is the ambipolar 
diffusion rate 
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with being the radar wavelength. Since the heating time scales are 
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so rapid (~10 ^ s) any measured change in y will be due to and not 
to T . or V . Hence the values of Y before and after heating can be 
compared to deduce the electron temperature increases, assuming that in 
the absence of heating the electron temperature is equal to a known 
neutral temperature. A similar comparison of values of the total Thomson 
scatter cross section would also reveal an electron thermal runaway, 
since “ (T +T.) The "after" measurement must be made within a few 

i. Si. 

milliseconds of turning off the heating pulse before the temperature 
rise decays. 

Let us now check our assumption that thermal conductivity is unim- 
portant. The appropriate scaling experiment is to have the radiowave 
heating rate much larger than the thermal conductivity rate y^> based 
on the beam radius QqZq. 
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In nondlmensional form, this equality becomes 
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In dimensional form it is 
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where = 3x10 s at an altitude of 100 km. Evidently, the neglect 
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of thermal conductivity is a good approximation everywhere in the E- 
'region,' ' 

Our analysis has shown that an appreciable fraction of the energy 
given to the electrons in the D— and E'-regions is given to vibrairional 
states and is resonantly transferred to CO^j which then radiates 4.3 U 
photons [Kumer and James, 1974; Kumer, 1977a, 1977b], Similarly,, energy 
given to 0^ vibrational energy is transferred to H^O, which radiates at 
6,5 [Houghton, 1967], Since vibrational excitation is limiting a 
thermal runaway in the D- and E-ragions, the efficiency for transforming 
absorbed radiowave energy into infrared photons must be high, probably 
exceeding: 80%. The- infrared radiation energy will be- split roughly 
equally between 4.3 ]i CO^ radiaion and 6.5 U H^O radiation. Therefore 
the amount of infrared radiation emitted can be estimated by computing 
the radiowave energy absorbed. 

The atteniiation length for the radiowave heating energy is defined- 
by 
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where the atmosphere scale height H = 10 km and we used v =10 s 


en 


because it is appropriate to the 1000°K electron temperature occurring 
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5 -3 

in a thermal runaway-. In the daytime E-region, n^ ~ 10 cm and the 
total attenuation of the radiowave energy is roughly 10% for V = 15 MHz. 
Obviously, with a lower frequency (5 MHz), the attenuation would almost 
complete. Hence, provided thermal runaway can be induced, roughly 
100 kW to 1 MW of power can be radiated as infrared photons. The 
overall infrared radiative power is roughly 10 W/cm -ster in the two 
bands having approximately 0.5 ]i bandwidth. These radiances are com- 
parable to the naturally occurring zenith radiance above 100 km and 
could be detected by a rocket-borne coherent detection instrument. 

Such an experiment would be an active way to test the theories of 
vibrational temperature and radiative transfer above 100 km. At night, 
the electron densities are roughly a factor of thirty lower, and the 
absorbed radiowave power is correspondingly less. 
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IONOSPHERIC MODEL 

The numerical model that is used to calculate the E~- and F-region 
ionospheric response to heating in the proposed satellite power station 
and Arecibo heating experiments is an extension of the model described 
by Roble [1975], Stolarski ^ [1975], and Roble [1976]. 

Figure 1 of Roble et al. [1977] is a block diagram showing the physical 
processes within the current model. A. list of the ion and neutral 
chemical reactions and the reaction- rates used within the model is also 
given in Roble et al. [1977] and Roble and Rees [1977]. Most of the 
physical and, chemical processes within the ionospheric model have been 
described in these earlier papers; therefore we will discuss here only 
the modifications and additions for this study.. However, since the 
heating experiments primarily affect the thermal structure of the iono- 
sphere, we will describe the processes involved in determining the 
ionospheric thermal balance in more detail - 

The time-dependent coupled electrcSh and ion energy equations 
within the model are solved simultaneously using the Newton iteration 
technique described by Hastings and Roble [1977] and Roble and Hastings 
[1977]. The electron heat sources include heating by fast photoelectrons 
in the daytime ionosphere, heat conduction from the magnetosphere, and 
resistive heating due to the absorption of radiowave energy beamed into 
the ionosphere as discussed in the second section of chis paper. Roble 
[1975] has shown that the thermal structure in the daytime ionosphere 
below about 250 km is maintained primarily by heating due to photoelectrons 
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and above that altitude by thermal conduction from the magnetosphere. 

The effects of thermal conduction in the electron gas become important 
above about 200 km, whereas below that altitude the electron tempera- 
ture is determined by a local balance of electron heating and cooling 
rates. The electron conductivity used in the analysis is given by 
Banks and Kockarts [1973] and Rees and Roble [1975], 

The ion heat sources are Coulomb coupling with electrons, chemical 
reactions [Stolarksi et al, , 1975], and heat conduction from the magneto- 
sphere. The ions lose energy in collisions with neutrals through polar- 
ization interactions and resonance charge exchange interactions. The 
ion-neutral loss rates are described by Banks and Kockarts [1973] and 
are listed in Table 3 of Rees and Roble [1975], Thermal conduction in 
the ion gas becomes important above about 300 km. 

The electron loss processes include energy loss in Coulomb colli- 
sions with ions as well as elastic and inelastic collisions with the 


neutrals. The elastic loss rates are given by (22.69) through (22,72) 
in Banks and Kockarts [1973]. The Inelastic loss rates include 
excitation of the fine structure levels of the ground state of atomic 
oxygen [(1) of Hoegy, 1976], excitation of atomic oxygen to the 0(^D) 
level [(61) of Rees and Roble, 1975], excitation of the vibrational 
levels of molecular nitrogen [(22,75) of Banks and Kockarts, 1973], 
and the excitation of the vibrational levels of molecular oxygen [(62) 
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Various processes within the model have been examined in detail in 
a recent comparison of model predictions with data obtained simultaneously 
by the Atmospheric Explorer-C satellite and the incoherent scatter radar 
station at Millstone Hill, Massachusetts [Roble et al, , 1977]. There was 
general agreement between model predictions and the satellite and inco- 
herent scatter radar measurements of electron and ion temperatures, 
electron and ion densities, 6300 1 airglow emission rates, NO densities, 
and photoelectron spectra in the daytime' ionosphere at the time of the 
satellite crossing over Millstone Hill. We use this same model to cal- 
culate the ionospheric response to radiowave heating from a satellite 
power station and from a Arecibo heating experiment, however, the simu- 
lation is for typical geomagnetic quiet conditions over an arbitrary 
midlatitude station. 

To simulate the effects on the ionosphere from the satellite power 
station and Arecibo heating experiment we consider two types of iono- 
spheric runs: steady state calculations in the midlatitude daytime 

ionosphere and diurnal calculations. In the daytime ionosphere the 
thermal and photochemical time constants are short, therefore the model 
is rtin to steady state [Roble et al. , 1977]. For these calculations, 
the neutral atmospheric temperature and composition profiles are speci- 
fied by the OGO-6 empirical model [Hedin et al. , 1974], the Hintereggar 
[1976] solar flux values are used, and typical daytime electron and ion 
heat fluxes and plasma flow from the magnetosphere are specified. These 
are the only values that are specified and the model is then run to 
steady state with all interactive processes being calculated self- 
consistently. The radiowave heating due to the satellite power station 
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and Arecibo heating facility are superimposed upon the normal daytime 
processes in this background ionosphere and the model is run to steady 
state. The differences between the two calculations give the. iono- 
spheric perturbations' due to radiowave heating. 

For the diurnal variation, a calculation similar to those- given by 
Roble [1975], Roble et al. [1976], and Roble et al. [1977] is performed 
for conditions typical in the geomagnetic quiet midlatitude ionosphere. 
The calculation is then repeated with the radiowave heating superimposed. 

The electron temperatures in the D- and lower E-regions of the iono- 
sphere are calculated on the basis of a local energy balance between 
electron heating' and -cooling- rates. The neutral atmospheric temperature 
profile and number density profiles of atmospheric constituents within 
this region are obtained from the U.S. Standard Atmosphere ' [NASA et al. , 
1976]. 
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IONOSPHERIC RESPONSE TO 5ABI0WAVE HEATING BY THE 
SATELLITE POWER STATION 


The satellite power station will beam 3 GHz, 10 GW radiowave energy 
from a geosynchronous orbit over the equator to a ground station some- 
where in the United States as shown in Figure 1. For our calculations, 
we arbitrarily consider two stations that collect the radiowave beam: 
Boulder, Colorado, and Houston, Texas. Although the two stations are 
geographically close, the radiowave beam intersects the geomagnetic 
field lines over them at different angles. As discussed in the second 
section, this results in a different' ionospheric heating rate over each ' 
station. The geomagnetic dip angle D is 70° for Boulder and 49“ for' 
Houston. From the geometry given in Figure 2 the angle S between the 
radiowave beam and the geomagnetic field line is 26,4° for Boulder and 
9° for Houston. The resulting values for the parameter given by (4), 
are 5.6 and 15.9 km for Boulder and Houston, repectively. 

The electron temperature response to the satellite power station 
radiowave heating calculated for the ionosphere over Boulder is shown in 
Figure 4. The ionospheric model calculates properties along individual 
geomagnetic field lines; therefore this figure is constructed using the 
results from several model runs. The peak radiowave heating occurs 
where the center of the beam intersects a geomagnetic field line. The 
radiowave heating distribution at this interaction is given by (3) and 
it is superimposed upon the normal daytime ionospheric heating rates. 

'The daytime ionospheric temperatures calculated outside of the microwave 
beam are indicated by the horizontal contour lines . 
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Below about 150 km, the electron temperature is determined by a 
balance between the local electron heating and cooling rates within the 
ionosphere. In the undisturbed daytime ionosphere the electron tempera- 
ture is nearly equal to the neutral gas temperature below about 120 km. 
The radiowave heating raises the electron temperature to a peak of about 
950°K throughout the D- and E-region of the ionosphere, causing a signi- 
ficant perturbation in the ambient ionospheric electron temperature 
along the radiowave beam. Above about 150 km, thermal conduction in the 
electron gas becomes important. Even though the bulk of the radiowave 
heating takes place at the beam’s intersection of the geomagnetic field 
line, thermal conduction in the electron gas redistributes this energy 
upward and it couples with the normal downward-directed heat flow from 
the magnetosphere in the daytime ionosphere. The two energy flows com- 
bine to increase the electron temperature along the geomagnetic field 
line. The maximum F-region temperatures occur along the geomagnetic 
field line that intersects the radiowave beam near the F^ peak. The 
region of enhanced electron temperature in the topside ionosphere follows 
the geomagnetic field line, whereas in the region below about 150 km it 
occurs along the radiowave beam. 

The calctilated electron temperature response in the ionosphere over 
Houston is shown in Figure 5, The angle between the radiowave beam and 
the geomagnetic field line over Houston is 9“ and, as a result, the F- 
region electron heating along a given geomagnetic field line that inter- 
sects the radiowave beam is greater over Houston than over Boulder, The 
calculated electron temperature is also greater over Houston than over 
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Boulder. However, the ionospheric area that is affected is smaller; 
over Houston- the radiowave- beam heats a narrow region more intensely 
Below about 150 lot, on the other hand, the calculated electron tempera- 
ture is similar in magnitude to the previous calculations, indicating 
local heating -and cooling along the radiowave beam throughout the D- 
and E-regions.- 

Maximum electron heating occurs when the radiowave beam is parallel 
to the geomagnetic field line.. The calculated electron temperature for 
this limiting case is shown in Figure 6, as are the calculated electron 
and ion temperature profiles for the undisturbed ionosphere. In the 
topside ionosphere the? calculated electron temperature is approximately 
4000®K along the geomagnetic field line that is coincident with the 
center of the radiowava beam. The electron temperature decreases toward 
both edges, of the radiowava beam. Below about 150 km,, the calculated 
electron temperature is similar to that in the previous two cases. 

The electron energy loss rates corresponding to- the undisturbed 
daytime ionosphere and to the case of maximum heating from the satellite 
power station (Figure 6), in which the radiowave beam is parallel to the 
geomagnetic field line, are shown in Figures 7a and 7b, respectively. 

In the undisturbed daytime ionosphere, the maximum- electron energy loss 
above about 200 km is due to Coulomb collisions with ions. Below about 
200 km, excitation of the vibrational levels of molecular nitrogen and 
molecular oxygen, the excitation of the ground levels of atomic oxygen, 
and the excitation of the rotational states of molecular nitrogen are 
all important processes' in electron energy loss. For comparison, the 
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total neutral gas heating rate is shown as the dashed curve in Figure 7a, 

In addition to the electron loss processes, the neutral gas heating con- 
sists of components due to ion chemistry [Stolarski et al. , 1975] and to 
heating from the photodissociation of molecular oxygen in the Schumann-Runge 
continuum [Roble, 1975], 

The electron loss rate corresponding to the case of maximum radio- 
wave heating when the satellite power station radiowave beam is parallel 
to the geomagnetic field line is shown in Figure 7b . Above about 300 km, 
the major electron loss process is Coulomb collisions with electrons and 
the excitation of atomic oxygen to the 0(^D) level. Below 300 km, the 
excitation of the molecular nitrogen vibrational levels becomes important. 
The rate of electron loss due to the excitation of the ground levels of 
atomic oxygen is most important between about 100 and 150 km. Finally, 
below 100 km the excitation of the vibrational and rotational levels of 
molecular nitrogen and molecular oxygen are the dominant electron energy 
loss processes. These results Indicate that the electron thermal runaway 
of the type discussed by Holway and Meltz [1973] and Meltz al. [1974] 
occurs below 100 km and that above that altitude thermal runaway can 
still occur but it is restrained by electron energy loss due to the 
excitation of the ground levels of atomic oxygen as discussed in the 
third section of this paper. The radiowave energy input to the ambient 
plasma in the ionospheric D- and E-regions is ultimately transferred to 
the neutral gas through collisions. 

Within the ionospheric F-region the radiowave heating produces 
enhanced airglow emission rates. Part of the atomic oxygen excited to 
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the 0C4)) level radiates as 6300 A. radiation and part is quenched by 
collision with molecular nitrogen and molecular oxygen. The calculated 
6300 1 volume emission rates in both the undisturbed daytime ionosphere 
and the maximum heating case are shown in Figure 8. The solid curve is 
the total 6300 k volume emission rate calculated in the daytime iono- 
sphere. The components include excitation of O(^) due to photodissocia- 
tion of molecular oxygen in the Schumann-Runge continuum, excitation due 
to photoelectron impact, excitation due to dissociative recombination 
of 0^, and thermal electron impact. The various excitation mechanisms, 
the excitation cross sections, chemical reaction rates, and quenching 
and' radiative coefficients used in the analysis have been discussed in 
detail by Roble et al. [1976]. The 6300 A volume emission rate peaks 
at 240 km and has a calculated dayglow intensity of 4.14 KR. The 
excitation component includes 1280 R due to Sch\imann-Runge photodissocia- 
tion, 1624 R due to photoelectron impact, 1191 R due to dissociative 
recombination, and 45 R due to thermal electron impact. The calculated 
6300 A volume emission rate for the maximum heating case is the dashed 
curve in Figure 8. The height of the 6300 A volume emission peak is 
near 260 km, with a total intensity of 13,1 KR. Almost all of the 
enhanced 6300 A airglow intensity in the daytime ionosphere is due to 
the impact of the hot ambient electrons on atomic oxygen along the 
radiowave beam. The dashed curves in Figure 8 give the calculated 
6300 A volume emission rate along various geomagnetic field lines for 
the Boulder case shown in Figure 4, The curves are for field lines 
where the radiowave beam intersects the geomagnetic field at 200, 250, 
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and 350 km. The calculated 6300 X intensities along these field lines 
are 4.41, 5.73, and 4.34 KR, respectively. 

The electron energy that excites molecular nitrogen vibrational 
states in the lower ionosphere is resonantly transferred to CO^, which 
then radiates 4,3 U photons [Kvimer and James, 1974; Kumer, 1977’a, 19.77b]. 
Similarly, energy given to molecular oxygen vibrational energy is trans- 
ferred to H^O, which radiates at 6.5 y. Hence, most of the radiowave 
energy that is transferred into the vibrational states of N 2 and 0^ is 
radiated to space as infrared photons. We have not done a detailed 
radiative transfer calculation to examine the magnitude of the enhanced 
4.3 y and 6.5 y radiations-; however, it is anticipated that the enhanced 
radiation is small compared to the black body spectrum radiated by the 
earth. Locally, however, the enhanced volxjme emissivity may be comparable 
to natural levels.. 

The electron density profiles calculated for the undisturbed daytime 
ionosphere and for the maximum heating case are shown in Figure 9. Below 
about 200 km, the electron density in the radiowave beam is enhanced by a 
factor of two, primarily because of the decreased reaction rate for dis- 
sociative recombination of the dominant molecular ions in the lower iono- 
sphere. The neutral and ion chemical reactions and reaction rates that 
are used in the' model are given in Table 1 of Roble et al. [1977]. The 
F 2 peak electron density decreases slightly; however, above the F 2 peak 
the electron density increases because of thermal expansion in the top- 
side ionosphere. Even though the electron temperature enhancement is 
large, the change in electron density in the F-region is relatively 
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small [Schiink and Walker, 1969]. The electron density profile was cal- 
culated three hours after the initiation of the radiowave heating and 
the response at the top of the ionosphere is time-dependent. However, 
it is unlikely that heating along any given field line will occur very 
long, since typical daytime ionospheric electric fields [Kirchhoff and 
Carpenter, 1976] will cause the geomagnetic field tube to drift out of 
the beam and be replaced' by another undisturbed geomagnetic - field tube. 

Within the high electron temperature region of the radiowave beam 
the vibrational levels of molecular nitrogen are excited by thermal 
electron impact [Newton et al. , 1974]. The increased vibrational tern- 
perature- of molecular nitrogen results, in an increased 0 + ^ NO -r N 

reaction rate [McFarland et al. , 1973], which' in turn causes a decrease 
in the F— region electron densities [Newton and Walker, 1975]. Our cal- 
culations have not included this process, which could result in lower F^ 
peak electron densities than those calculated in Figure 9, 

The satellite power station will beam radiowave energy into the 
ionosphere during the day and night. To determine the diurnal iono- 
spheric response to this heating, the diurnal ionospheric properties 
are calculated for a case in which there is no radiowave heating and 
compared to a case in which the radiowave heating is superimposed. We 
consider the maximum heating case when the radiowave beam is parallel 
to the geomagnetic field liner The diurnal variations in electron 
temperature and density, calculated over a typical midlatituda iono- 
spheric station, are shown in Figures 10a and 11a, respectively. The 
electron density increases rapidly after sunrise and reaches a maximian 
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C 

of 10 cm at 260 km near 1400 LT. ^ The electron density in the F- 
region decreases thereafter,' reaching minimum values just before sun- 
rise. Below the F 2 peak the electron density decreases to nighttime 
values rapidly after sunset-. The bulge in the electron density near 
1300 LT is due to an enhanced downward flow of o"*" from the magneto- 
sphere that is arbitrarily assumed in the calculation. 

The diurnal variation of the calculated electron temperature shows 
an inverse relationship with electron density in the daytime ionosphere. 
The enhanced electron temperatures in the topside ionosphere in the 
morning and evening hours are due to the enhanced heat flow rates from 
the magnetosphere that are assumed in 'the calculation. The diurnal 
variation of the ionosphere is typical of midlatitude ionospheric condi- 
tions. 'We have not attempted to model a specific day, as described by 
Roble [1975], Roble et al. [1976], and Roble et al. [1977], but merely 
calculate some typical ionospheric variation from which the affects of 
the radiowave heating by the satellite power station can be evaluated. 

As discussed in the second section, the radiowave heating of the 
ambient plasma is proportional to the electron density and electron- 
neutral collision frequency. The calculated diurnal variation of the 
radiowave heating for the case when the radiowave beam is parallel to 
the geomagnetic field line is shown in Figure 12. The radiowave heating 
of the ambient plasma has a diurnal variation that is similar to the 
diurnal variation of electron density. Maximum heating occurs in the 
lower D- and E-regions where the electron-neutral collision frequency is 
higho At F-region heights maximum heating occurs in the daytime iono- 


sphere 
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The calciilated diurnal variations of electron temperature and 
density when the satellite power- station radiowave beam is parallel to 
the geomagnetic field line over the ionospheric station are shown in 
Figures 10b and lib, respectively ». The calculated electron temperature 
is considerably enhanced from typical ionospheric, conditions, and below 
about 300 km the diurnal variation is small. The electron temperature 
variation reflects not only variations in the electron heating rate due 
•to radiowave heating but also variations in the topside heat flow rate 
in the electron gas that is assumed from the magnetosphere. The calcu- 
lated electron temperature due to the radiowave heating is relatively 
constant throughout the day because the electron heating rates and elec- 
tron loss rates both depend upon the electron density and electron 
neutral collision frequency in the same manner. 

The calculated diurnal variation in electron density for' the case 

when the radiowave beam is heating the ionosphere is shown in Figure lib. 

The electron density at the F^ peak is slightly lower than in the case 

of no radiowave heating; however, in the topside ionosphere the plasma 

density has increased. This is primarily due to thermal expansion in • 

the topside ionosphere in response to the enhanced electron temperatures. 

The high electron temperatures also give rise to enhanced 6300 A airglow 

% 

emission rates throughout the day. 
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IONOSPHERIC RESPONSE TO RADIOWAVE HEATING BY THE 
PROPOSED ARECIBO HEATING EXPERIMENT 

As discussed in the second section of this paper, the ionospheric 
heating rates of a 3 GHz, 10 GW radiowave beam from a satellite power 
station can be closely simulated by a 15 MHz, 3 MW heating experiment 
from a ground-based station. In this section, we present calculations 
of the ionospheric response to various heating experiments from an 
assumed ground-based station at Arecibo. The ionospheric calculations 
are similar to those described in the previous section; however, the 
geometry is modified to simulate a heating experiment from a ground- 
based station. The assumed beam shape, power flux, and ionospheric 
volume heating rate associated with a 15 MHz, 3 MW heating experiment 
have been described in the second section. 

The calculated ionospheric electron temperature response to a 
15 MHz, 3 MW heating experiment is shown in Figure 13. Below about 
150 km, the ground-based radiowave beam heats the electrons locally to 
high temperatures. In the D-region the ambient electron temperature 
within the beam is raised by an order of magnitude from ambient back- 
ground electron temperatures . The calculated electron temperatures in 
the E-region are slightly lower than D-region temperatures, and above 
about 150 km, thermal conduction along the geomagnetic field line 
causes the major heating response to occur in the direction along the 
field line rather than along the beam. Since the volume heating rate 
is proportional to the electron density and electron- ion and electron- 
neutral collision frequencies, the peak heating within the F-region 
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occurs near the peak and decreases with altitude in the topside iono- 
sphere. Within the F— region the calculated electron temperatures for 
the ionospheric heating experiment are similar to chose calculated for 
the satellite power station. In the D-region, however, the ground-based 
heating experiment causes a considerably larger perturbation in the 
ambient electron temperature than the satellite power station radiowave 
beam. The large perturbation could be used to probe D-region processes, 
as is proposed in the discussion section. 

A similar calculation is .made for a 15 MHz , 1 MW heating experi- 
ment, and the resiilts are shown in Figure 14. The calciilated tempera- 
ture response is smaller than in the previous case; however, the tem- 
perature pattern is similar. For comparison with the satellite power 
station ionospheric effects, the 15 MHz, 3 tlW heating experiment produces 
similar effects at F-region altitudes whereas the 15 MHz, 1 MW heating 
experiment produces siiailar effects in the D— and E-regions, The cal- 
culated electron temperature for the maximum heating case, when the 
15 MHz, 3 MW radiowave beam on Arecibo is directed along the geomagnetic 
field line, is shown in Figure 6. The calculated F-region temperature 
is similar to the maximxim heating produced by the satellite power 
station; however, the D- and E-region ionospheric temperatures are 
greater. The calculated electron energy loss rates for this maximum 
heating experiment are shown in Figure 7c. In the F-region the electron 
energy loss rates are similar to the loss rates calculated for the 
satellite power station.. In the D- and E-regions of the ionosphere the 
major energy loss processes are the excitation of the vibrational states 
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of molecular oxygen and molecular nitrogen followed by the excitation 
of the rotational levels of the same molecules. This is just the 
opposite of the loss rate sequence for the satellite power station 
because the vibrational excitation loss rates increase rapidly with 
temperature and the 15 MHz, 3 MW heating experiment produces greater 
electron temperatures in the D- and E-regions than the satellite power 
station. The increased vibrational excitation loss rates will also 


give rise to increased 4,3 y and 6.5 y radiation because of the effi- 
cient exchange of vibrational energy between molecular nitrogen and 
carbon dioxide and molecular oxygen and water vapor in the region below 
about 100 km., The high electron temperatures also excite the 0,^(a^A ) 

and 0 (b^Z ) levels by thermal electron impact, thus giving rise to 
^ S 

enhanced airglow emissions at 1.27 y and 8644 1 , respectively [Noxon, 
1970, 1975]. In the F-region the quenching of 0(^D) by O 2 results in 
an enhanced 0_ (b^Z ) emission of about 100 R which could be detected 
from the ground at night. The 15 MHz, 3 MW Arecibo heating experiment 
transfers enough radiowave energy to the ambient plasma to double the 
neutral gas heating rate within the radiowave beam from its value under 
typical daytime conditions. 


The calculated electron density response to heating by the 15 MHz, 


3 MW Arecibo heating experiment is shown in Figure 15. Below about 

150 km, the molecular ions O 2 and NO are the dominant ions. The O 2 

ion dissociative recombination rate used in the calculation varies as 

for T < 1000°K and for T > 1000°K, The NO’*’ ion dissocia- 

e e ■ e e 

tive recombination rate varies as T for T > 1000®K and for 

e e e 
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< lOGO^IC. Both reaction rates are given in Table 1 of Roble et al, 
[1977], along with other important temperature-dependent rate coefficients 
involved, in the O 2 and. NO chemistry,. The increased electron temperature 
in the D— and E-region results in an increase in the eleccron density 
that is 2.5 times greater in the vicinity of the radiowave beam thaii in 
the region outside the beam. There is a slight decrease in the electron 
density in the vicinity of the peak', with an increase above the peak. 
This effect is due to thermal expansion in the topside ionosphere and 
it is time-dependent. The calculations shown in Figure 15 are for con- 
ditions three hours after the radiowave heating was imposed. It is 
unlikely that the field tube wotild remain longer in the radiowave beam 
because the ambient electric field should replace field tubes within 
the beam at a more rapid rate. The region of enhanced electron density 
in the topside ionosphere is along the geomagnetic field line in the 
same’ direction as the region of enhanced electron temperature. 

The Arecibo radiowave heating experiment produces the largest elec- 
tron temperature perturbations in the ionospheric D- and E-regions . The 
electron temperature in these regions is a strong function of the radio- 
wave heating power as shown in Figure 16. In the D-region, the electron 
temperature increases almost linearly with the heating power until about 
2.5 Mff, when the losses due to vibrational excitation of molecular oxygen 
and molecular nitrogen become dominant and change the slope of the elec- 
tron temperature curve. In the D-region the atomic oxygen number density 
is small so that the major- loss rates are due to rotational and vibra- 
tional excitation of molecular oxygen and nitrogen. Above about 90 km. 
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the loss due to the excitation of the fine structure levels of atomic 
oxygen also become important. This condition is shown in Figure 16, 
where it is seen that the slope of the electron temperature curves with 
respect to radiowave heating power is less for altitudes above 90 km. 
than for altitudes below that. 
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DISCUSSION 


The previoiis sections have shown that the satellite power station 
radiowave beam perturbs the ionospheric properties from normal back- 
ground values and that these effects can be simulated by an Arecibo 
heating facility. A 15 MHz, 3 MW Arecibo heating facility produces 
effects at F— region heights similar to those produced by a satellite 
power station. The heating and ionospheric response for both experi- 
ments depend upon the angle between the radiowave beam and the geo- 
magnetic field line. Within the E-region the ionospheric temperature 
response is. about two times greater for the Arecibo heating experiment 
than for the satellite power station. The effects are comparable at 
D- and E-region heights if the power of the Arecibo heating experiment 
is reduced to 1 MW. Thermal runaway occurs in the D-region because the 
main loss processes at low thermal energies is the electron impact 
excitation of rotational states of and 0^. The thermal runaway, 
however, is limited at high teaperatures by the electron impact excita- 
tion of the vibrational levels of and 0^. Within the E-region electron 
cooling due to the excitation of the fine structure of atomic oxygen also 
limits thermal runaway. It should be possible to use the measured D- 
and E-region electron densities and temperatures to infer the atomic 
oxygen density in the^ altitude region above about 90 km, 

Airglow measurements, should also prove to be a diagnostic measure 
of aeronomic processes in the high temperature region. Enhancements in 

the volume emission. rate of 0„ (a^A ) and 0„(b^S ) radiation at 1,27 ]i and 

^ a ^ S 
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.8644 1. could be detected in addition to the enhanced 6300 A emission , at 
F-region heights. The infrared' volume emission rates at 4.3 y due to 
CO 2 emission and .at 6.5 y due to H 2 O emissions within the radiowave 
beam could be detected locally from space. Measurements of the doppler 
width of the 6300 -A airglow emission- line with a Fabry-Perot --inf erero- 
meter can be used to determine the neutral gas temperature at F-region 
heights. Similarly measurements of the rotation bands of 0.^(b^S ) can 
be used to determine the neutral gas temperature in the D-region and the 
doppler width of the atomic oxygen green line gives the neutral gas 
temperature near 97 km.. 

Finally, the D- and E-region electron density and temperature varia- 


tions can be used to derive information on reaction rates of various 


aeronomic processes, such as the recombination rates. 

Our calculations show plasma density enhancements along the geo- 
magnetic field line in the topside ionosphere over the ionospheric 
region heated by the radiowave beam. These enhancements would presum- 
ably extend along the geomagnetic field tube toward the conjugate hemi- 
sphere and perhaps interfere with natural processes- occurring within the 
plasmasphere. 

Evidently, the most dramatic changes caused by a satellite power 

station occur in the ionospheric D- and E-regions where thermal runaway 

occurs. To induce the same effect at Areciho, a heater with a power 

4 2 

aperture product PA > 7x10 MW m is required. (We have .assumed 
A^_^ = 0.5A.) Therefore an underdense heating experiment at a frequency 
V = 15 MHz, 2 i-lW of radiated power, and efficient use of the Arecibo 
dish should' simulate the ionospheric heating aspects of a satellite 
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power station. But we should clearly point out that any facility which 
falls short of these specifications will fail to produce the thermal 
runaway and hence will not simulate the satellite power station. And, 

I 

because the electron temperature will remain low, no interesting 
molecular airglow lines will be excited that can be used for analyzing 
D- and E-region dynamical, chemical, and aeronomical effects. 
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Fig. 1 


.Fig. 2 


FIGURE CAPTIONS 


a) Geometry of the satellite power station beam with respect 
to. the earth’s magnetic field B which has a dip angle D. 
The angle 9 is the latitude and B is the angle between the 
microwave beam and the earth’s magnetic field, 

b) Diagram determining the altitude extent of the heating 
along a magnetic field line ^?hich passes through the center 
of the microwave beam. 

Geometry of an underdense radiowave heating experiment. The 
altitude interval A over which the radiowave beam intersects .a 
given magnetic field line is small in comparison with the size 
of the ionosphere itself. The solid curve represents a pro- 
jection of the heating rate per unit volume along the magnetic 
field line onto the vertical direction, and AT is a scher- 
matic representation of the temperature change along the 
magnetic field line . Specifically, lines A and B connect a 
particular temperature change AT^ with the spatial location 
where that temperature change occurred. Note that thermal 
conductivity produces temperature changes well outside the 
radiowave heating beam. Also sketched is the altitude resolu- 
tion Ah appropriate to a photometer line of sight. 

Relation between the critical values of A [see (30), where 
it is . denoted by A^] and a [see. (31)] at which thermal runaway 

occurs. The formula for A is -(ZS),, 

c 


Fig. 3. 
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Fig. 4. Contours of calculated electron temperature (°K) over Boulder 

for radiowave heating by the satellite power station. The long 
dashed lines give the radiowave beam direction and the light 
solid lines indicate the geomagnetic field direction. 

Fig. 5. Same as Figure 4, except for heating over Houston. 

Fig. 6. Vertical profiles of electron temperature (°K) calculated for 
the case where the radiowave beam is parallel to the geomagne- 
tic field line giving maximum electron heating. The solid 
curves give the calculated electron and ion temperatures in 
the undisturbed daytime ionosphere. The dashed curve is the 
calcula'ted electroh temperature due to maximum heating by the 
satellite power station, and the long-and-short-dashed curve is 
the electron temperature due to maximum heating by the 15 MHz, 

3 MW Arecibo heating experiment. The daytime neutral gas tern- 
perature profile is labeled T^ and the dashed T^ curve is the 
calculated ion temperature for both the satellite power 
station and the Arecibo heating experiment. 

Fig. 7. Calculated vertical profiles of the total neutral gas heating 

-3 -1 

rate (dashed curve, total-n), eVcm s , the total electron 
loss rate in collisions with neutrals (solid curve, total-e), 
and various electron loss components. The cuirve labels indi- 


cate electron loss rates; N 2 J,, 0^^, and 0^ are due to elastic 


collisions with N^, 0^, and 0, respectively; and are 
due to excitation of rotational structure of N 2 and O 2 ; N 2 ^ 


and 02 y are due to excitation of vibrational levels of O 2 and 
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?xg. 8. 


Fig. 


9. 


Fig, 10. 


N 2 ; 0(^P) and O(^) are due to excitation of atomic oxygen 

fine structure and atomic oxygen to the ^ level. 0^(b^S ) 

1 11 
and 0^(a A ) are due to excitation of 0„ to the Z and A 
2 g 2 g g 

levels by electron impact, and e— i indicates electron-ion 
Coulomb collisions.- The cases are (a) loss rates in undis- 
turbed daytime ionosphere, (b) maximum heating by the satellite 
power station, and (c) maximum heating by the 15 MHz, 3 MW 
Arecibo heating experiment.. 

The calculated vertical profile of the 6300 A voliane emission 
—3 —1 

rate (photons cm s ) . The solid curve is the volume emis- 
sion rate in' the undisturbed daytime ionosphere, and the long- 
and-short-dashed curve is the emission rate calculated for 
maximum heating by the satellite power station. The dashed 
curves are volume emission rates calculated over Boulder along 
geomagnetic field lines that intersect the radiowave beam at 


200, 250, and 350 km, respectively. 

—3 

Calculated vertical profiles of electron density (cSi ) . The 
solid curve is the electron density calculated in the undis- 
turbed daytime ionosphere and the dashed curve is the electron 
density curve determined three hours after the' maximum radiowave 
heating due to the satellite power station is turned on. 

Contours of the diurnal variation of electron temperature (°K) 
calculated (a) in the undisturbed daytime ionosphere, and 
(b) with the maximum heating rata of the satellite power 
station superimposed. 
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Fig. 11. 


Fig. 12. 


Fig. 13. 


Fig. 14. 
Fig. 15. 
Fig. 16. 


-3 

Contours of the diurnal variation of electron density (cm ) 
calculated (a) in the undisturbed daytime ionosphere, and 
(b) with the maximum heating rate of the satellite power 
station superimposed. 

-3 -1 

Diurnal variation of electron heating rate [log^^CQ^, ergs cm s )] 
in the case of maximum radiowave heating by the satellite power 
station. 

Contours of calculated daytime electron temperature (°K) over 
Arecibo in response to the 15 MHz, 3 ionospheric heating 
experiment. The solid lines with arrows are geomagnetic field 
lines. The heating experiment beam is Indicated by the diverg- 
ing solid lines. 

Same as Figure 13, except for a 15 MHz, 1 MW heating 
experiment, 

-3 

Same as Figure 13, except electron density (cm ) for 
the 15 MHz, 3 MW heating experiment. 

Relationship between electron temperature and heating power 
for an Arecibo heating experiment at various altitudes in the 
D— and E-regions. 
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APPENDIX 2 

IONOSPHERE-MICROWAVE INTERACTION SUBGROUP 
0. Garriott, W. E. Gordon, L, Duncan 

The purpose of this subgroup was to discuss the possi- 
bility for scaled laboratory studies of ionosphere-microwave 
interactions-, in particular thermal self-focusing instabilities. 
The JSC Chamber A facility is well-suited to the large scale 
sizes of these interactions and we need to deteirmine if the 
other plasma parameters can be scaled appropriately. 

The transmitter system was . proposed to be a microwave 
horn producing a beam incident on the laboratory plasma. The 
fact that this is not a plane wave, as assumed in the theory, , 
needs to be considered. To avoid a standing wave configuration 
in the plasma, and thus to avoid plasma density striations un- 
related to self- focusing, the microwave radiation should be 
absorbed after passing across the chamber. The question arose 
as to whether absorbing material exists which can be put into 
a vacuum environment. 

In a later conversation it was suggested that, under 
the proper geometry, microwave reflection could be tolerated. 
However, if microwave absorption can be achieved in a vacuum, 
it is the better solution because it suppresses the unwanted 
standing waves. 

We understand that a plasma generator could produce a 
large plasma, as required for Chamber A. Other methods, such 
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as using an electron beam, .could produce plasma instabilities 
much stronger than the self- focusing mechanism. 

The question of plasma diagnostics was not addressed. 

In the general discussion group, Langmuir probes were the 
most commonly mentioned density diagnostic. Light and radio 
wave scattering are other possibilities. The usefulness of 
each method is dete3onined by its sensitivity, which then must 
be related to the An/n of the density striations . The time 
scales and saturation levels of these density fluctuations 
has not been discussed. 

Although close attention . has been given to the perpen- 
dicular scale size of the density striations / we have 

not considered the effects of finite parallel scale sizes, 
limited by the laboratory dimensions . 

It was suggested that we might try parametric excita- 
tion in an overdense plasma (a situation much closer to pre- 
vious work in the ionosphere) as a means of building confidence 
in our laboratory technique . 

The major problem in producing thermal self-focusing 
instabilities in a plasma laboratory is the scaling of the 
plasma parameters . 

2 

We assume a threshold power of ~ 1 kw/m as a 

reasonable laboratory microwave transmitter output. A beam 
diameter of 3-5 meters is in keeping with the chamber dimen- 
sions. So as to produce several striations within the plasma, 
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we require X < 2m (striation width = X /2) . To avoid 

X X 

2 2 

parametric interactions , ^ 5 ^ critical ' 

operating frequency and ^critical plasma critical 

frequency. 

using - (6 X cos^e p'l 

yields f^rit ~ ^^crit ^ ^ 10^^] , neglecting 9 term; 


f .075 GHz [n . . ~ 7 x 10 ] , 8 = 89®; 

crxt crxt 

f . . ~ 1.8 GHz [n . . 4 X 10^*^] , 9 = 65° ~ vertical 

incidence for 
Chamber A. 


10 3 

We understand a plasma density of 10 electrons/cm is difficult 
to obtain . 

A question was raised in the general group discussion con- 
cerning the T^ dependence of X^. This dependence, working through 
the electron collision frequency, does not appear in the expres- 
sion given above because of an ionospheric- like plasma approxima- 
tion. With explicit T^ dependence included. 



(3.44 X cos^e ±1 


The plasma density and T^ are related through the electron mean 
free path, which is limited by the condition 

X ^ « chamber dimensions, 

mrp 

For reasonable values of T , this expression reduces to the 

0 .to 

previous equation. 
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One possibility now being expiored is allowing neutrals 
to. act as the- primary souce of electron collisions, replacing 
electron-ion collisions. Then, 


X" 


(3 X 10^^)cos^e f“^. (T /n^)p“?; 

o crit e' n th 


To maintain low conductivity across field lines, we require 
« f^^ (collision frequency « gyro frequency) . 


For 


T 


j. 

9 


5 f 

crxt 
10^ °K 


2 m 
65° 




then n ~ 4.1 x 10 
n 


13 




V- ~ 2.7 X 10^ (f ~ 10^)' 
en ce 


Notice that f can be increased by using an applied field. 

This capability also was suggested by the beam-plasma subgroup. 
It is probably advantageous to maintain an ionospheric-like 
condition that f ./f ““4. This approach appears promising, 

]p0 C0 ■“ 

but additional thinking is required. 

In order to determine if Chamber A is well suited for 
use in studying thermal self- focusing instabilities, we need 
to answer the following questions : 
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1 ) Transmitter 

(a) What are the amplitude variations in the near field? 

(b) Can absorbing material be used in a vacu\m? 

(c) Are there^ hardware problems v/ith the microwave system 
in a vacuum? 

2 ) Plasma production 

What n , n , T , L (scale height) are readily produced? 
110 6 

3) Plasma diagnostics 

(a) What are the operating characteristics of a Langmuir 
probe? 

(b) Can we use microwave radio wave, or laser scattering? 

(c) What other diagnostics might be effective? 

4) Theoretical studies 

(a) What is the saturation prediction for An/n ? 

(b) What is the effect of electron-neutral collisions? 

(c) What is the effect of finite j| dimensions to the 
striations ? 

5) What are the amplitude and shape of applied B-fields? 

6) What is necessary for parametric excitation? (i.e., repeat 
questions 1-5) 
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